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ABSTRACT 

We present the Evolution of molecular Gas in Normal Galaxies (EGNoG) survey, an observational 
study of molecular gas in 31 star- forming galaxies from z = 0.05 to z = 0.5, with stellar masses of 
(4 — 30) X 10^^ M0 and star formation rates of 4 — 100 M© yr~^. This survey probes a relatively 
un-observed redshift range in which the molecular gas content of galaxies is expected to have evolved 
significantly. To trace the molecular gas in the EGNoG galaxies, we observe the CO (J = 1^0) and 
CO (J = 3 -^ 2) rotational lines using the Combined Array for Research in Millimeter- wave Astronomy 
(CARMA). We detect 24 of 31 galaxies and present resolved maps of 10 galaxies in the lower redshift 
portion of the survey. We use a bimodal prescription for the CO to molecular gas conversion factor, 
based on specific star formation rate, and compare the EGNoG galaxies to a large sample of galaxies 
assembled from the literature. We find an average molecular gas depletion time of 0.76 ± 0.54 Gyr 
for normal galaxies and 0.06 =b 0.04 Gyr for starburst galaxies. We calculate an average molecular 
gas fraction of 7-20% at the intermediate redshifts probed by the EGNoG survey. By expressing the 
molecular gas fraction in terms of the specific star formation rate and molecular gas depletion time 
(using typical values), we also calculate the expected evolution of the molecular gas fraction with 
redshift. The predicted behavior agrees well with the significant evolution observed from z ^ 2.5 to 
today. 



1. INTRODUCTION 

In the past decade, molecular gas observations have 
begun probing the high redshift universe in a sys- 
tematic way using increasingly powerful millimeter in- 
struments. The picture that is emerging at redshifts 
1 — 2 is similar in some respects to what we see in 
the local universe. Sub-mm galaxies (SMGs) are ob- 
served to be undergoing extreme starbursts as a result 
of major mergers (e.g. Tacconi et al. 2008; Engel et al. 
2010), equivalent to local ultra-luminous infrared galax- 
ies (ULIRGs) (e.g. Sanders et al. 1986; Solomon et al. 
1997; Downes & Solomon 1998). Normal star- forming 
galaxies (SFGs) at high redshifts, akin to local spi- 
rals, are becoming accessible as well. Recent works 
(Tacconi et al. 2010; Daddi et al. 2010) suggest that z ~ 
1 — 2 star-forming galaxies (with star formation rates 
(SFRs) of ~ 50 — 200 M© yr~^) are scaled up versions 
of local spirals, forming stars in a steady mode (not 
triggered by interaction), despite hosting star formation 
rates at the level of typical local starburst systems like 
luminous infrared galaxies (LIRGs) and ULIRGS. 

While galaxies classified as LIRGs or ULIRGs (by their 
infrared luminosities only) have typically been associated 
with starbursting and merging systems by analogy to 
galaxies in the local universe, it is becoming clear that 
this connection does not hold at high redshifts. Morpho- 
logical studies find that while 50% of local LIRGs show 
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evidence of major mergers (Wang et al. 2006), that frac- 
tion appears to decrease toward high redshifts: Bell et al. 
(2005) find that more than half of intensely star-forming 
galaxies at z ^ 0.7 have spiral morphologies and fewer 
than 30% show evidence of strong interaction. Further, 
the typical star formation rate of SFGs increases toward 
higher redshift. SFGs have been observed to obey a tight 
relation between stellar mass and star formation rate 
out to 2; ~ 4 (e.g. Brinchmann et al. 2004; Noeske et al. 
2007; Elbaz et al. 2007; Daddi et al. 2007; Pannella et al. 
2009; Daddi et al. 2009; Magdis et al. 2010). This 'main 
sequence' of SFGs evolves with redshift, shifting to higher 
star formation rates at higher redshifts. 

The increase in the star formation rate is mirrored 
in the molecular gas fraction (/mgas = ^mgas/(^* + 
^mgas), where M^gas is the molecular gas mass (in- 
cluding He) and M* is the stellar mass) of these sys- 
tems. While studies of local spirals (e.g. FCRAO sur- 
vey. Young et al. 1995; BIMA SONG, Regan et al. 2001, 
Heifer et al. 2003; HERACLES, Leroy et al. 2009) find 
average molecular gas fractions of ~5%, observations of 
high redshift SFGs suggest molecular gas fractions of 20- 
80%, an order of magnitude higher than local spirals. 
Unfortunately, the past 8 Gyr of the evolution of /mgas 
remains relatively unprobed, with only one study of nor- 
mal SFGs between z = 1 and z = 0.05 (Geach et al. 
2011, at z = 0.4). 

To fill in this observational gap in redshift, we have 
carried out the Evolution of molecular Gas in Normal 
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Redshift 


Redshift 


Sample 


Parent 


M* Range 


^FR Range 


Observed 


^obs 


OARMA 


Bin 


Range 


Size 


Sample 


(Mo) 


(Mo yr-l) 


Transition 


(GHz) 


Configuration 


A 


0.05 - 0.10 


13 


SDSS 


(4.0- 16)xl0^^ 


3.4 - 41 


CO( J = 1^0) 


107 


C 


B 


0.16 - 0.20 


10 


SDSS 


(6.3- 32) X 10^0 


47 - 106 


CO( J = 1^0) 


98 


D, some C 


C 


0.28 - 0.32 


4 


SDSS 


(16- 32) X 10^0 


39 - 65 


CO( J = 1^0) 
CO( J = 3^2) 


88 
266 


D 
D, E 


D 


0.47 - 0.53 


4 


zCOSMOS 


(4.0- 5.5)xlOio 


62 - 74 


CO( J = 3^2) 


230 


D, E 



Table 1 

Summary of the EGNoG survey. 



Galaxies (EGNoG) survey-^, a key project at the Com- 
bined Array for Research in Mihimeter-wave Astronomy 
(CARMA).^ The EGNoG survey uses rotational transi- 
tions of the carbon monoxide (CO) molecule to trace the 
molecular gas content of 31 star- forming galaxies from 
z = 0.05 to 0.5 (the CO(J =1^0) line for galaxies 
at z = — 0.3 and the CO(J = 3^2) line for galax- 
ies at z = 0.3 — 0.5). The survey includes four galaxies 
at z ~ 0.3 (the gas excitation sample) observed in both 
the CO(J =1^0) and CO(J = 3^2) lines, more 
than doubling the number of normal SFGs at z > 0.1 in 
which CO line ratios have been measured. The results 
from the gas excitation sample are presented separately, 
in Bauermeister et al. (2013). 

In this paper, we present the data for the entire EG- 
NoG sample and use these data to constrain the evolution 
of the molecular gas fraction in galaxies at intermediate 
redshift. The paper is organized as follows. In Section 2 
we describe the selection of the EGNoG survey galaxies 
and the properties of the sample. In Section 3 we describe 
the observations and data reduction (see Appedix B for a 
detailed description of the data reduction and Appendix 
C for polarization measurements of calibrators 0854+201 
and 1058+015). We present the CO maps, total CO lu- 
minosities, and derived molecular gas masses in Section 4 
and discuss the morphologies of the low-redshift portion 
of the sample as well as the non-detections of the high- 
redshift portion. In Section 5, we present a compilation 
of data from the literature and analyze the behavior of 
the molecular gas depletion time and molecular gas frac- 
tion in starburst and normal galaxies at low and high 
redshift. We provide some concluding remarks in Sec- 
tion 6. 

Throughout this work, we use a ACDM cosmology with 
(h, Qm, ^a) = (0.7, 0.3, 0.7). We use the following terms 
for star- forming galaxies at low and high redshift: 

• Star-forming galaxy refers to any galaxy which 
is actively forming stars. Quantitatively, these 
galaxies lie on or above the main sequence of star- 
forming galaxies. This term encompasses both nor- 
mal star- forming and starburst galaxies. 

• Normal star-forming galaxy or SFG refers to 
star-forming galaxies that lie on the main sequence 
(within some scatter). While the acronym SFG 
stands for star- forming galaxy, we use it (in keep- 
ing with other authors) to refer to main sequence 
galaxies only, which comprise the bulk of the star- 
forming galaxy population. 

^ More information and data products of the EGNoG survey can 
be found at http://carma.astro.umd.edu/wiki/index.php/EGNoG 

■^ CARMA is a 3-band, 23-element millimeter interferometer 
jointly operated by the California Institute of Technology, Uni- 
versity of California Berkeley, University of Chicago, University of 
Illinois at Urbana-Champaign, and University of Maryland. 



• Starburst galaxy refers to a galaxy undergoing 
a period of enhanced star formation. We give a 
quantitative definition in Section 2.1. 

• LIRG, ULIRG (luminous Infrared galaxy, ultra- 
luminous Infrared galaxy) denotes a galaxy with a 
total infrared luminosity between 10^^ and 10^^ Lq 
(LIRG) or > 10^2 L^ (ULIRG). Local LIRGs and 
ULIRGs tend to be starbursts, but at high-redshift, 
normal star-forming galaxies would be classified as 
LIRGs or ULIRGs based on their infrared luminosi- 
ties alone. 

• SMG (sub-mm galaxy) refers to the population of 
galaxies at z > 1 which have 850-/im fluxes > 4 
mJy. SMGs tend to be starbursts (often as a result 
of a merger, akin to local ULIRGs). 

2. THE EGNOG SAMPLE 

Utilizing the CO(J = 1^0) and CO(J = 3 ^ 2) in 
the 3 mm and 1 mm bands of CARMA, it is possible to 
achieve continuous redshift coverage from z = to z ^ 
0.5. The EGNoG survey observed 31 galaxies spanning 
the accessible redshift range in 4 redshift bins, A-D. The 
redshift range, sample size, parent sample and range of 
M^ and SFR are given for each redshift bin in Table 1. 

Sample galaxies in redshift bins A-C are drawn from 
the main spectroscopic sample of the Sloan Digital 
Sky Survey (SDSS), Data Release 7 (York et al. 2000; 
Strauss et al. 2002; Abazajian et al. 2009). Spectro- 
scopic redshifts are from David Schlegel's spZbest files 
produced by the Princeton- ID code, specBS^. The stel- 
lar masses and SFRs of galaxies in the SDSS DR7 
are provided by the Max-Planck-Institute for Astro- 
physics - John Hopkins University (MPA-JHU) group 
(http : //www . mpa-garching . mpg . de/SDSS). Stellar 

masses are derived by fitting SDSS ugriz photometry to 
a grid of models spanning a wide range of star forma- 
tion histories. This method is found to compare quite 
well with the Kauffmann et al. (2003) methodology us- 
ing spectral features (more detail on this comparison is 
found on the website above). Star formation rates are de- 
rived by fitting the fluxes of no less than 5 emission lines 
using the method described in Brinchmann et al. (2004). 
Both stellar masses and star formation rates are derived 
using a Bayesian analysis, producing probability distri- 
butions of each quantity for each galaxy. The probability 
distributions of the stellar mass and SFR for galaxies in 
redshift bins A-C are given in Figure Al in Appendix 
A. We take the median of the distribution, with errors 
indicated by the 16th and 84th percentile points (ilcr 
for a Gaussian distribution). In cases where a dupli- 
cate SDSS source exists (as a result of SDSS automated 
source-finding), we take the average of the two median 

^ See http://spectro.princeton.edu/ for more information 
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EGNoG 

Name 


SDSS 
Identification 


RA 


Dec 


z 


log(M*/M0) 


SFR 

(M© yr-i) 


sSFR 

(Gyr-i) 


Al« 
A2 
A3 
A4 
A5 
A6 
A7 
A8 
A9 
AlO 
All 
A12 
A13 


SDSS J234311. 26+000524.3 
SDSS J231332.46+133845.3 
SDSS J233455.23+141731.0 
SDSS J085504.16+525248.4 
SDSS J085307.26+121900.8 
SDSS J084549.66+573239.3 
SDSS J211527.81-081234.4 
SDSS J135751. 77+140527.3 
SDSS J105733.59+195154.2 
SDSS J141601. 21+183434.1 
SDSS J100559.89+110919.6 
SDSS J111150. 65+281147.7 
SDSS J221938.11+134213.9 


23:43:11.257 
23:13:32.468 
23:34:55.239 
08:55:04.171 
08:53:07.258 
08:45:49.662 
21:15:27.817 
13:57:51.775 
10:57:33.589 
14:16:01.216 
10:05:59.890 
11:11:50.654 
22:19:38.115 


+00:05:24.319 
+13:38:45.301 
+14:17:31.093 
+52:52:48.320 
+12:19:00.827 
+57:32:39.329 
-08:12:34.441 
+14:05:27.317 
+19:51:54.274 
+18:34:34.171 
+11:09:19.688 
+28:11:47.767 
+13:42:13.890 


0.096941 + 0.000007 
0.081196 + 0.000008 
0.062168 + 0.000010 
0.089592 + 0.000028 
0.080978 + 0.000011 
0.085259 + 0.000005 
0.090469 + 0.000006 
0.099188 + 0.000027 
0.077338 + 0.000007 
0.055121 + 0.000008 
0.076055 + 0.000006 
0.098171+0.000012 
0.083503 + 0.000012 


10.74+°o:J°9 
10.99l°:}°o 

10.92l°o:fo 
ILOllg;}? 
10.801°;°} 

10.86lg:J°9 

10.55±g:S°9 
10.75lg:g^9 
10.72±g:g^8 
10.55lg:S^ 

ii.06i°:}J 


28.4+ti 
21.5t^:° 
49.9lt;^ 
42.3tej 

40.9112,1 

3.7111 

5.11?:! 
s.gl^^-o 

3.41?:^ 

,, 7+12.8 
11. '-6. 8 


0.54 
0.29 
0.20 
0.60 
0.42 
0.56 
0.57 
0.05 
0.10 
0.09 
0.17 
0.10 
0.10 


Bl 
B2 
B3 

B4« 
B5 
B6 
B7 
B8 
B9 

BIO 


SDSS J223528.63+135812.6 
SDSS J002353.97+155947.8 
SDSS J100518.63+052544.2 
SDSS J105527.18+064015.0 
SDSS J115744.35+120750.8 
SDSS J124252.54+130944.2 
SDSS J091426.24+102409.6 
SDSS J114649. 18+243647.7 
SDSS J134322.28+181114.1 
SDSS J130529.30+222019.8 


22:35:28.638 
00:23:53.973 
10:05:18.640 
10:55:27.188 
11:57:44.348 
12:42:52.548 
09:14:26.239 
11:46:49.182 
13:43:22.288 
13:05:29.304 


+13:58:12.619 

+15:59:47.757 
+05:25:44.225 
+06:40:15.025 
+12:07:50.792 
+13:09:44.215 
+10:24:09.649 
+24:36:47.703 
+18:11:14.135 
+22:20:19.867 


0.183193 + 0.000017 
0.191773 + 0.000020 
0.165737 + 0.000012 
0.173078 + 0.000013 
0.182661 + 0.000010 
0.174661+0.000013 
0.176180 + 0.000018 
0.176579 + 0.000021 
0.178004 + 0.000019 
0.190227 + 0.000010 


11.43+°:}^ 
ii.29ig:}? 

10 79+0-09 

iU.^y_Q Qg 

11 n2+o-ii 
10.83l°:?l 

-I -I .P' + O.IO 

11 nr+o-ii 

-,-, of:+0.12 
ii.oo_Q -,^2 
10 qfi+0-09 


oo 9+129.9 

54.61^5:^ 

47 1+55.9 
^'..^-23.3 
40 7+82.8 
^y. '-28.7 
CI 1+19.3 
01. 1-14.5 
cq ,+ 15.9 

6i.5l|?:^ 

loe.ilif,^ 

+84.8 

75.011^,^ 


0.33 

0.28 
0.76 
0.47 
0.29 
0.87 
0.22 
0.89 
0.30 
0.82 


Cl^ 
C2 
C3 
C4 


SDSS J092831. 94+252313.9 
SDSS J090636.69+162807.1 
SDSS J132047.13+160643.7 
SDSS J133849.18+403331.7 


09:28:31.941 
09:06:36.694 
13:20:47.139 
13:38:49.189 


+25:23:13.925 
+16:28:07.136 
+16:06:43.720 
+40:33:31.748 


0.283020 + 0.000022 
0.300622 + 0.000010 
0.312361+0.000014 
0.285380 + 0.000015 


11.241°:^° 
11.201°:?! 
11.461°:?^ 
ii.26i°:i« 


38.7lfJ 

57 r + 90.1 

64.9i^fo' 
50.5lf,-i 


0.22 
0.37 
0.23 
0.28 


Dl^ 
D2^ 
D3^ 
D4^ 


SDSS J100055.81+015703.8 
SDSS J100052.41+014833.0 
SDSS J095939.07+022249.6 
SDSS J095900.61+022833.0 


10:00:55.82 
10:00:52.41 
09:59:39.07 
09:59:00.62 


+01:57:03.84 
+01:48:32.92 
+02:22:49.94 
+02:28:33.20 


0.474100 + 0.000367 
0.527900 + 0.000367 
0.470800 + 0.000367 
0.477800 + 0.000367 


io.6oi°:i^ 

10.73l°:i^ 
10.74l°:i^ 


62.2l-i 

78.7i^^-3 
74.4t^?-^ 


1.55 
1.47 
1.47 
1.36 



Table 2 

Basic information for EGNoG galaxies. 
^ Indicates duplicate source in SDSS. The average value is reported for 2;, M* and SFR. 
^ Bin D sources are selected from COSMOS. COSMOS identifications for galaxies D1-D4 are zCOSMOS 811469, 

COSMOS 2019408 and zCOSMOS 840823, respectively. 



^COSMOS 811543, 



values and use the lowest (highest) 16th(84th) percentile 
value to indicate the negative (positive) error. 

The higher redshift portion of our sample (redshift bin 
D) is drawn from the Cosmic Evolution Survey (COS- 
MOS; 2 square degrees at RA ^ 10 h, Dec '^ 2°; 
Scoville et al. 2007) which has imaging and photomet- 
ric redshifts for all the galaxies in the survey. Spectro- 
scopic redshift s are available for many sources (zCOS- 
MOS; Lilly etal. 2009). We did not use the SDSS 
for this redshift bin due to the poor coverage of star- 
forming galaxies at this redshift. Lilly et al. (2009) re- 
port an average accuracy of 110 km s~^ for the spectro- 
scopic redshifts, independent of redshift. Stellar masses 
come from spectral energy distribution (SED) fitting by 
Bundy et al. (2010) and star formation rates are based 
on SED fits provided by Ilbert et al. (2010). Typical er- 
rors are 0.15 dex for the stellar masses and 0.3-0.4 dex 
for the star formation rates. 

The redshifts, stellar masses and star formation rates 
for the entire EGNoG sample are given in Table 2. 



2.1. Sample Selection 

The EGNoG galaxies were selected from the parent 
samples (SDSS and COSMOS) to be as representative 
as possible of the 'main sequence' (MS) of star-forming 
galaxies. The main sequence is the tight correlation be- 
tween M* and SFR that has been observed over a large 
range of redshifts: e.g. z ^ (Brinchmann et al. 2004), 
^ - 0.2 - 1 (Noeske et al. 2007), z - 1 - 2 (Elbaz et al. 
2007; Daddi et al. 2007; Pannella et al. 2009), ^ - 3 - 4 
(Daddi et al. 2009; Magdis et al. 2010) (see also the sum- 
mary of recent results in Dutton et al. 2010). The sSFR 
increases with redshift, tracing the increase in the cosmo- 
logical dark matter (and cold baryons, fuel for star for- 
mation) accretion rate onto halos (Bouche et al. 2010). 
Building on these observations, a few authors have at- 
tempted to describe the main sequence relation at all 
redshifts with one form. Figure 1 shows the func- 
tional forms of the specific star formation rate (sSFR= 
SFR/M*) as a function of redshift reported by three 
papers: Bouche et al. (2010), Karim et al. (2011) and 
Elbaz et al. (2011). Bouche et al. (2010) suggest that 



Bauermeister et al. 



J..U 


' ' ' ' . 
















• ' .••'L ^ 




• ' >**^ 




. * " •■' ^ 








•'■'•:;'"^'*'^ ^^ 


0.5 


> '.^ 




,' •> .^^^y^^y^ 




* -^ ^^^^^^''^ 




•^ ^^^^^o^"^ 




-^ .^^^.t^^ .y^ 




- .^ .^^i^^^^y^ 




* ^ .^^<y^ .^^ 




'' ^' .^r^y^^y^ ^ 




- <^- .-^^X^^^^-- — 


0.0 


' v^ ^^^^^^^^^^^'^^ 








- ^y ^^y .j^^ 




.* ^y 4^^/^ ^^^ 




♦ # .' ^V^ >^^ 




♦ V .•' V>^ .y^ 


0.5 


ySy^ 








Elbaz2011 






-- Bouche2010 logM,=10 




w/ 


Bouche2010 logM,=ll 




-■- Karim2011 logM, =9. 8-10.2 


1.0 


4^ • 


Karim2011 logM, =10.2-10. 6 " 






- - Karim2011 logM, =10. 6-11.0 




^ X 


Karim2011 logM,>11.0 

^ this work logM, =11 


1.5 




- - this work logM, = 10 






1 



0.0 



0.5 



1.0 



1.5 



2.0 



2.5 



Figure 1. Specific star formation rate (sSFR) versus red- 
shift for different stellar masses. Three published forms are 
shown (Elbaz et al. 2011 in blue, Bouche et al. 2010 in red, and 
Karim et al. 2011 in green) as well as the relation adopted for this 
work in black. 

existing literature data are well-fit by the form sSFR 
(X M-^-^{\ + zf-^ for z = - 2. Karim et al. (2011) use 
a deep 3.6 /im selected sample of galaxies from the COS- 
MOS field, with average SFRs determined from stacked 
1.4 GHz radio continuum emission. Like Bouche et al. 
(2010), they find the sSFR to be a function of stel- 
lar mass, with galaxies with lower stellar masses hav- 
ing higher specific star formation rates. Specifically, the 
authors find sSFR ex M~^-^(\ + zf-^ (on average) for 
star- forming galaxies at z = — 3 (in Figure 1, we plot 
the sSFR fit for each of the Karim et al. (2011) mass bins 
individually). On the other hand, Elbaz et al. (2011) ar- 
gue that observations at all redshifts are consistent with a 
sSFR independent of stellar mass. The authors use far- 
infrared observations from Herschel Space Observatory 
in conjunction with existing data on the GOODS-North 
and GOODS-South fields to fit the sSFR versus redshift, 
finding sSFR (Gyr-^) ex t~^^l^,^ for z = - 3 (where 
^cosmic is the time since the Big Bang in Gyr). The sSFR 
forms from Bouche et al. (2010) (red) and Karim et al. 
(2011) (green) are shown for different stellar masses in 
Figure 1. As Elbaz et al. (2011) argues that the sSFR is 
independent of sSFR, only one curve (blue) is plotted. 

For this work, we are only concerned with the behav- 
ior of sSFR(2;) out to z ?^ 2.5, but note that there is 
evidence for a fiattening of this relation at z > 2 — 2.5 
(e.g. Stark et al. 2009; Gonzalez et al. 2010; Reddy et al. 
2012). For z < 2.5, we adopt a sSFR relation which 
roughly agrees with the three forms described above. 
Specifically, we define 



sSFRMs(Gyr-^) = 0.07(1 + z) 



3.2 



M, 



1011 Mc7 



-0.2 



(1) 



In Figure 1, this is plotted in black for M* = 10^^ Mq 
(dashed) and lO^^ Mq (solid). We use this representative 
form of the sSFR to describe the main sequence of star- 



forming galaxies throughout this work. 

To differentiate 'starburst' (SB) galaxies from nor- 
mal main sequence galaxies, we adopt the quantitative 
definition given by Rodighiero et al. (2011), a study of 
PACS/Herschel observations of COSMOS and GOODS- 
South galaxies at z = 1.5 — 2.5. Rodighiero et al. (2011) 
fit the main sequence population with a Gaussian dis- 
tribution and define starbursts as outliers of this distri- 
bution. They find deviations from the Gaussian main 
sequence starting at a sSFR four times that of the main 
sequence. Therefore, we define starbursts as having 



sSFRsB > 4 X sSFRms 



(2) 



In the selection of the EGNoG sample, we apply the 
following criteria in each redshift bin in order to iden- 
tify non-interacting, star-forming galaxies lying as close 
to the main sequence as possible. At z < 0.3 (bins A- 
C), star- forming galaxies were selected (rejecting sources 
with AGN) using the cut from Kauffmann et al. (2003) 
in the BPT line-ratio diagram (Baldwin et al. 1981). 
We did not extend this criterion to bin D at z ~ 0.5 
since the Kauffmann et al. (2003) dataset only covers 
z = 0.02 — 0.3. Obviously interacting galaxies were ex- 
cluded via visual inspection of the SDSS or COSMOS 
optical images. However, some interacting galaxies may 
remain in the sample due to the difficulty of identification 
at the modest resolution of the SDSS images. 

Practical considerations imposed the following further 
constraints. From the SDSS dataset, we only used spec- 
troscopically targeted galaxies, for which spectroscopic 
redshifts, stellar masses and SFRs were available from 
the MPA-JHU group. We required a spectroscopic red- 
shift so that the error in the redshift is small enough to 
ensure that CO emission would be captured within the 
observed bandwidth. We excluded galaxies with SFRs 
below a minimum value estimated from the instrument 
sensitivity in each redshift range, assuming a molecu- 



lar gas depletion time (rrr 



M^gas/SFR) of 2 Gyr 



(Leroy et al. 2008), a Milky Way-like CO conversion fac- 
tor (e.g. Dame et al. 2001) and (for bins C and D) a 
CO(J = 3 ^ 2) to CO(J =1^0) line ratio (rsi) 
of 0.6 (Aravena et al. 2010). Note that this is close to 
rsi = 0.5, found in the EGNoG bin C galaxies at z ^ 0.3 
(see Bauermeister et al. 2013). In bin A, we used two 
SFR cutoffs since this bin is composed of two subsam- 
ples: high-SFR galaxies (EGNoG A1-A6) observed as a 
pilot study for the EGNoG survey, and a lower- SFR cut- 
off sample (EGNoG A7-A13) that is more representative 
of the main sequence. Finally, in order to restrict our 
selection to the main sequence at the SFRs probed, we 
required a minimum stellar mass of lO^^"^ M©. 

In each redshift bin, the EGNoG galaxies were selected 
randomly from the subset of the parent sample that sat- 
isfied the constraints discussed above. Figure 2 shows 
stellar mass versus SFR in each of the four redshift bins. 
EGNoG galaxies are indicated by the red points with 
error bars. The vertical and horizontal dotted lines indi- 
cate the minimum M* and SFR, respectively, required 
for sample selection in each redshift bin. Galaxy CI 
lies below the minimum SFR for bin C because it is a 
dupHcate source in the SDSS catalog and was selected 
based on the higher SFR. The average SFR of the dupli- 
cate entries lies below the minimum SFR for bin C. The 
blue shading indicates the (logarithm of the) density of 
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Figure 2. Stellar mass versus SFR in each redshift bin. The redshift bin is indicated in the upper left corner of each panel. Red points 
with error bars show EGNoG galaxies. Dotted vertical and horizontal lines indicate the minimum M* and SFR, respectively, required for 
sample selection in each redshift bin (see text for full description). Blue shading shows the logarithm of the density (in the M*-SFR plane) 
of star-forming galaxies from the parent dataset in the specified redshift range (given in the lower right corner). The redshift ranges plotted 
in bins B-D are slightly larger than the EGNoG redshift ranges, so that more data points may be included to better capture the behavior 
of the main sequence at each redshift. The solid black line indicates the main sequence (Equation 1) at the average redshift of the bin. 
The starburst criterion (Equation 2) is indicated by the black dashed line. 



points in the M* - SFR plane for all star-forming galax- 
ies (spectroscopic targets only for the SDSS data, bins 
A-C) in the redshift range indicated in the lower right of 
each panel. The redshift ranges plotted for bins B-D are 
slightly larger than the EGNoG redshift ranges so that 
more points may be included to better capture the be- 
havior of the main sequence at each redshift. The blue 
shading corresponds to approximately 100,000 galaxies 
at z = 0.05 - 0.1, 15,000 galaxies at z = 0.15 - 0.21, 
1300 galaxies at z = 0.25 — 0.35, and 14,000 galaxies 
at z = 0.43 — 0.57. The main sequence of star-forming 
galaxies at each redshift is indicated by the solid black 
line, with the starburst cutoff indicated by the dashed 
black line. From bin A to bin C, the low- mass, low- SFR 
end of the main sequence becomes sparsely sampled as 
the number of spectroscopically targeted, star-forming 
galaxies available in the SDSS decreases. This effect 
is most dramatic in the bin C panel, where the sample 
of spectroscopically-targeted star- forming galaxies in the 
SDSS is an order of magnitude smaller than in the bin B 
redshift range. Conversely, the small area covered by the 
COSMOS survey results in a sparsely sampled high-mass 
end of the main sequence in the bin D panel. 

The result of our sample selection is that the EGNoG 
galaxies generally lie at the high-M*, high-SFR end of 
the main sequence of star- forming galaxies. Our star- 
burst criterion suggests that roughly half of the galaxies 
in each of bins A and B, and all of the bin D galaxies 
are starbursts. All four bin C galaxies lie below the star- 
burst cut. In each redshift bin, the EGNoG galaxies lie 
roughly within the width of the scatter around the main 
sequence. The starburst classification of some of the EG- 
NoG galaxies will be considered further in Section 5. 

3. CARMA OBSERVATIONS 

Each of the 31 EGNoG galaxies was observed in at least 
one rotational transition of the CO molecule: galaxies 
in bins A, B and C were observed in CO (J = 1^0) 
(z/rest = 115.3 GHz, z/obs ^ 107, 98 and 88 GHz re- 
spectively); galaxies in bins C and D were observed in 
CO(J = 3^2) (z/rest = 345.9 GHz, u^hs ^ 266 and 
230 GHz respectively). At each frequency, each galaxy 
was observed over several different days. Each dataset 
includes observations of a nearby quasar for phase cali- 
bration (taken every 15-20 minutes), a bright quasar for 
passband calibration and either a planet (Uranus, Nep- 



tune or Mars) or MWC349 for flux calibration (in most 
cases). 

The reduction of all observations for this survey was 
carried out within the EGN^ data reduction infras- 
tructure (based on the MIS pipeline; Pound & Teuben 
2012) using the Multichannel Image Reconstruction, Im- 
age Analysis and Display (MIRIAD, Sault et al. 2011) 
package for radio interferometer data reduction. Our 
data analysis also used the mi ri ad-python software 
package (Williams et al. 2012). The data were flagged, 
passband-calibrated and phase calibrated in the stan- 
dard way. Final images were created using invert with 
options=inosaic to properly handle and correct for the 
three different primary beam patterns. All observations 
are single-pointing. We describe the CO (J = 1^0) and 
CO (J = 3^2) observations individually below. A full 
description of the data reduction and flux measurement 
is given in Appendix B. 

3.1. CO(J = 1^0) 

The CO (J =1^0) transition lies in the 3 mm band 
of CARMA (single-polarization, linearly polarized feeds) 
for galaxies in bins A, B and C. The line was observed 
with at least three overlapping 500 MHz bands, cover- 
ing ^ 4200, 4600 and 5000 km s"^ total, at 35, 39 and 
42 km s~^ resolution for galaxies in bins A, B and C, 
respectively. 

Bin A galaxies were observed during three time pe- 
riods: October to November 2010, April to May 2011 
and February 2012. Ah data were taken in CARMA's C 
configuration, with 26 — 370 m baselines yielding a typ- 
ical synthesized beam of 2.0'' x 1.5'' at 107 GHz. These 
observations are sensitive to spatial scales up to ^ 22", 
which is sufficient for most galaxies in the sample but 
may resolve out some large scale structure in the largest 
galaxies. However, since molecular gas is centrally con- 
centrated, we expect any under-estimation of the flux to 
be minimal. (In D and E configuration in the 3mm band, 
we are sensitive to sufficiently large spatial scales so that 
we expect to recover all of the flux.) Each galaxy was ob- 
served for 2-3 hours (time on-source), yielding final data 
cubes with a typical rms noise of 5 — 10 mJy beam~^ in 
35 km s~^ channels. 

^ http : //carma . astro . -iimd . edu/wiki/index . php/EGN 
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Bin B sources were observed from August to November 
2011 and April 2012 in CARMA's D configuration, with 
11 — 150 m basefines yielding a typical synthesized beam 
of 4.9'' X 3.9'' at 98 GHz. Supplementary observations 
were made of Bl, B2, B3 and B7 in CARMA's C array in 
February 2012, which resulted in a synthesized beam of 
3.0" X 2.3" in the final maps combining data from both 
array configurations. Each galaxy was observed in D 
configuration for approximately 7 hours (time on-source), 
resulting in a typical rms noise of 2.5 mJy beam~^ in 39 
km s~^ channels in the final cube. 

Bin C observations were made from August to Novem- 
ber 2011 in CARMA's D configuration, with a typical 
synthesized beam of 4.8" x 3.9" at 88 GHz. Each galaxy 
was observed for 20 to 30 hours (time on-source), yielding 
final cubes with a typical rms noise of 1.2 mJy beam~^ 
in 42 km s~^ channels. 

The fiux scale in each dataset is set by the fiux of 
the phase calibrator, which is determined from the fiux 
calibrator. For the fiuxes used, see Appendix B. 

Of the 27 galaxies observed in the CO(J =1^0) 
line, we detected 24 easily, with 3 non-detections in bin 
A. The ixv-spectra of the detected galaxies are shown in 
Appendix D. 

3.2. C0(J = ?>^2) 

The CO(J = 3^2) transition lies in the CARMA 

1 mm band (dual-polarization, circularly polarized feeds) 
for galaxies in bins C and D. We again observed the line 
with at least three overlapping 500 MHz bands, covering 
?^ 1500 and 1700 km s"^ total, at 14 and 15 km s"^ res- 
olution, for galaxies in bins C and D respectively. 

Bin C observations were carried out in CARMA's E 
configuration during August 2011 and D configuration 
during April 2012. Source C4 was observed entirely in 
E configuration, while the other three sources were ob- 
served mostly in D configuration. The E configuration 
has 8— 66 m baselines yielding a typical synthesized beam 
of 3.2" X 2.5" at 266 GHz. The D configuration has 
11 — 150 m baselines yielding a typical synthesized beam 
of 1.7" X 1.5" at 266 GHz. Each galaxy was observed for 

2 to 6.5 hours (time on-source), yielding final cubes with 
a rms noise of 5-10 mJy beam~^ in 42 km s~^ channels. 

Bin D observations were carried out in CARMA's E 
and D configuration during August 2011 and April to 
June 2012, respectively. In E configuration, a typical 
synthesized beam is 4.2" x 3.6" at 230 GHz. In D con- 
figuration, a typical synthesized beam is 2.5" x 1.7" at 
230 GHz. Galaxies Dl and D2 were observed in both 
D and E configuration for approximately 10 hours (time 
on-source) total, yielding final cubes with a rms noise of 
~ 7 mJy beam~^ in 48 km s~^ channels. Galaxies D3 and 
D4 were observed in D configuration for approximately 

3 hours (time on-source), yielding final cubes with a rms 
noise of ~ 10 mJy beam~^ in 48 km s~^ channels. 

In bin C, sources C2, C3 and C4 were detected at the 
^ 5cr level in the C0( J = 3^2) line. However, CI, with 
its wide velocity profile (observed in the CO (J = 1 -^ 
0) line), was only marginally detected and we provide 
an upper limit on the CO (J = 3 -^ 2) fiux. The uv- 
spectra for the bin C galaxies are shown in Appendix D. 
In bin D, we do not significantly detect any of the four 
galaxies, and list only upper limits. The bin D galaxies 



are discussed further in Section 4.3 

4. RESULTS 

4.1. CO Emission Maps and Molecular Gas Masses 

Table 3 presents the quantities we derive from the CO 
observations. The CO line luminosity is calculated from 
the line fiux {Sqo in Jy km s~^, calculated as described 
in Section B.2) following 



L'co = 3.25 X 10^ Sco vl 



.(1 + ^)" 



(3) 



(see the review by Solomon & Vanden Bout 2005), 
where i^obs is in GHz and rcom is the comoving distance 
in Mpc. The units of L^q are K km s~^ pc^. We report 
the measurement error for Sqq- The error reported for 
Lqo includes both the measurement error of ^co and a 
30% systematic error, added in quadrature (see Section 
B.2 for more details). 

The center velocity, '^^center, is the fiux- weighted aver- 
age velocity of the galaxy-integrated spectrum ('u = 
at the redshift in Table 2). The error reported is the 
standard deviation of the '^center values found with the 
three fiux measurement methods and different channel 
averaging described in Section B.2. The reported veloc- 
ity width (AF) is the full width of the emission, where 
'source' velocity channels are selected by eye. We give 
the velocity width of a single channel as the error. 

In the case of a non-detection (galaxies A4, A6, A12, 
Dl, D2, D3, and D4; indicated by ^ in Table 3), we 
estimate a 3cr upper limit on ^co from the noise in the 
channel maps and a maximum velocity width motivated 
by detections at that redshift: 600 km s~^ in bin A and 
400 km s-^ in bin D. For the CO(J = 3^2) line in 
galaxy CI, while the CO (J = 3 -^ 2) channel maps did 
not show evidence of a source upon visual inspection, 
an integrated spectrum made of a circular region 4.5" in 
radius at the center of the image suggests a 3cr detection. 
We calculate an upper limit on the line fiux (and related 
quantities) from this spectrum, over the velocities of the 
CO(J =1^0) line in CI. 

The molecular gas mass (M^gas) is calculated from 
L'qq according to 



M 



mgas 



1.36 aco ^c 



CO(J-(J-l)) 



/ m 



(4) 



where aco is the CO (J = 1^0) luminosity to H2 
mass conversion factor (akin to Xco) and the factor of 
1.36 accounts for Helium to yield the total molecular 
gas mass, rji is the ratio of the CO line luminosities: 
Tj\ = L'qq{J -^ J — 1)/L'qq(1 ^0). In the calculation 
of the molecular gas mass from the CO (J = 3 -^ 2) lumi- 
nosity (bin C and D observations), we use r^i =0.5, ob- 
served in the EGNoG bin C galaxies (Bauermeister et al. 
2013). We use an aco value selected according to the 
starburst classification described by Equations 1 and 2. 
For galaxies classified as starburst, we use the value ob- 
served in local ULIRGs: aco(ULIRG) = 0.8 M© (K 
km s~^ pc^)~^ (Scoville et al. 1997; Downes & Solomon 
1998). For all other galaxies, we use a Milky Way-like 
conversion factor, aco (Milky Way) = 3.2 M© (K km s~^ 
pc^)~^ (e.g. Dame et al. 2001). This classification (SB 
or norm), and the molecular gas mass calculated using 
the corresponding aco is reported for each galaxy in Ta- 
ble 3. There is of course considerable uncertainty in the 
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Name 


J up 


'S'co 


'^center 


Ay 




^CO 




sSFR 


-/Wmgas 


/mgas 


Tdep 






(Jy km s-i) 


(km s~ 


^) 


(km s~ 


^) 


(109 K km s 


-' Pc2) 


class 


(lO^Mo) 




(Gyr) 


Al 




10.75 ± 2.99 


33.7 ± 


27.3 


428.0 ± 


71.3 


4.76 ± 


1.95 


SB 


5.18 ± 


2.12 


omtVot 


n -17+0.09 
^•-L'^-0.08 


A2 




14.40 ± 2.55 


30.5 ± 


12.6 


281.2 ± 


35.2 


4.44 ± 


1.55 


norm 


19.32 ± 


6.73 


0.171°:°^ 


PI ^0+0.27 


A3 




17.94 ± 2.07 


22.8 ± 


6.8 


552.5 ± 


34.5 


3.22 ± 


1.03 


norm 


14.01 ± 


4.50 


o.i2tS;g^ 


omtr^l 


A4« 




5.98 ± 1.98 


0.0 ± 


0.0 


600.0 ± 


0.0 


2.25 ± 


1.01 


SB 


2.45 ± 


1.09 


0.031°:°} 


0.05l°:°o'. 


A5 




23.22 ± 3.90 


-60.3 ± 


9.8 


492.0 ± 


35.1 


7.12 ± 


2.45 


SB 


7.75 ± 


2.66 


r. p,7 + 0.03 


0-18i°o:o7 


A6« 




3.65 ± 1.23 


0.0 ± 


0.0 


600.0 ± 


0.0 


1.24 ± 


0.56 


SB 


1.35 ± 


0.61 


0.021°:°} 


n 04+0-02 

^•^^-0.02 


A7 




6.72 ± 1.52 


49.2 ± 


31.9 


319.0 ± 


106.3 


2.58 ± 


0.97 


SB 


2.81 ± 


1.05 


omtTol 


p p,7+0.03 


A8 




4.52 ± 1.35 


-61.0 ± 


20.5 


357.3 ± 


71.5 


2.10 ± 


0.89 


norm 


9.14 ± 


3.87 


o.o5i°:°3 


-1 99+0.70 
-L-^^_0.70 


A9 




5.68 ± 1.42 


-7.1 ± 


17.7 


420.3 ± 


70.1 


1.59 ± 


0.62 


norm 


6.92 ± 


2.70 


r. -1^ + 0.06 


-1 0^+0.91 


AlO 




7.88 ± 3.16 


4.3 ± 


17.0 


274.4 ± 


68.6 


1.11 ± 


0.56 


norm 


4.83 ± 


2.42 


0-08-0.04 


n Q4+0.54 


All 




6.89 ± 0.98 


-10.5 ± 


11.2 


419.8 ± 


70.0 


1.86 ± 


0.62 


norm 


8.09 ± 


2.69 


O.ut'o'ol 


o.9it°:i 


A12^ 




3.50 ± 1.16 


0.0 ± 


0.0 


600.0 ± 


0.0 


1.59 ± 


0.71 


norm 


6.92 ± 


3.09 


O.wt'ofr 


2.0311:^5 


A13 




21.89 ± 4.12 


-27.3 ± 


8.9 


563.5 ± 


35.2 


7.15 ± 


2.53 


norm 


31.12 ± 


11.02 


no-|+0.08 


9 «7+l-63 
^•'^'^-1.63 


Bl 




11.35 ± 0.62 


-23.1 ± 


4.0 


615.4 ± 


38.5 


18.60 ± 


5.67 


norm 


80.95 ± 24.68 


r. r,o+0.08 

U.ZO_Q Q7 


n Q9+0-60 


B2 




7.32 ± 0.40 


25.5 ± 


11.7 


581.1 ± 


38.7 


13.20 ± 4.03 


norm 


57.45 ± 


17.52 


n, 90+0. 08 

U.ZO_Q Q7 


1 0^+0-54 


B3 




4.49 ± 0.35 


-23.1 ± 


6.6 


378.9 ± 


37.9 


5.97 ± 


1.85 


SB 


6.50 ± 


2.01 




n 14+0.08 


B4 




4.81 ± 0.20 


-25.6 ± 


5.2 


457.6 ± 


38.1 


6.99 ± 


2.12 


SB 


7.61 ± 


2.30 


n 07+0-03 


n if;+o-26 


B5 




6.61 ± 0.22 


16.1 ± 


2.5 


307.6 ± 


38.5 


10.70 ± 


3.23 


norm 


46.57 ± 


14.06 


o.2itg:g^ 


P) q-|+0.38 
'^•^-'^-0.38 


B6 




3.52 ± 0.31 


-133.6 ± 


16.4 


343.7 ± 


38.2 


5.22 ± 


1.63 


SB 


5.68 ± 


1.78 


o.o8t°:S| 


o.iot°:gl 


B7 




5.97 ± 0.30 


-40.9 ± 


4.9 


497.0 ± 


38.2 


9.01 ± 


2.74 


norm 


39.21 ± 


11.93 


n 19+0-04 


0.641°:- 


B8 




2.25 ± 0.33 


-34.8 ± 


20.4 


382.5 ± 


76.5 


3.41 ± 


1.14 


SB 


3.71 ± 


1.24 


n ns+o-oi 


0.031°:°^ 


B9 




3.91 ± 0.27 


-18.6 ± 


5.9 


459.5 ± 


38.3 


6.03 ± 


1.86 


norm 


26.24 ± 


8.08 


n 11+0-04 


0.39l°:f, 


BIO 




5.25 ± 0.22 


27.1 ± 


2.6 


270.9 ± 


38.7 


9.28 ± 


2.81 


SB 


10.10 ± 


3.06 


n 10+003 


n -.0 + 0.05 
U'-L'^-O.OS 


CI 




2.05 ± 0.34 


-84.9 ± 


18.1 


542.2 ± 


41.7 


8.25 ± 


2.83 


norm 


35.90 ± 


12.31 


n -17+0.06 


^ qo+2.09 


C2 




2.35 ± 0.10 


-6.7 ± 


1.2 


253.7 ± 


42.3 


10.70 ± 


3.24 


norm 


46.57 ± 


14.11 


0.23l°:J? 


p. 0-1+0.39 
0-81-0.39 


C3 




4.39 ± 0.09 


-24.7 ± 


5.7 


384.0 ± 


42.3 


21.70 ± 


6.53 


norm 


94.44 ± 28.40 


0.25l°:« 


1.451°:- 


C4 




3.01 ± 0.11 


17.3 ± 


13.0 


292.5 ± 


41.8 


12.30 ± 


3.72 


norm 


53.53 ± 


16.18 


p r,o + 0.10 


i.06i°:« 


CI" 


3 


9.02 ± 2.65 


-49.9 ± 


7.3 


542.3 ± 


23.4 


4.03 ± 


1.69 


norm 


35.08 ± 


14.73 


p -17+0.07 
^•-'-'^-0.07 


0.9ltl°l 


C2 


3 


8.39 ± 1.12 


19.5 ± 


8.2 


169.2 ± 


42.3 


4.25 ± 


1.40 


norm 


36.99 ± 


12.15 


o.i9i°:j° 


U.o4_Q 32 


C3 


3 


13.46 ± 1.56 


-14.0 ± 


7.8 


426.7 ± 


42.7 


7.39 ± 


2.38 


norm 


64.32 ± 20.69 


o.i8tS;S^ 


0.99l°:g 


C4 


3 


11.88 ± 1.22 


33.5 ± 


5.0 


250.7 ± 


41.8 


5.41 ± 


1.72 


norm 


47.09 ± 


14.93 


^21+0-10 


0.93l°:« 


Dl" 


3 


3.36 ± 1.12 


0.0 ± 


0.0 


400.0 ± 


0.0 


4.38 ± 


1.96 


SB 


9.53 ± 


4.27 


n 1Q+0-09 


0.151°:- 


D2" 


3 


2.73 ± 0.91 


0.0 ± 


0.0 


400.0 ± 


0.0 


4.45 ± 


2.00 


SB 


9.68 ± 


4.34 


n -17+0.08 
^•-•-'-0.07 


n 1S+0-18 


D3" 


3 


3.39 ± 1.13 


0.0 ± 


0.0 


400.0 ± 


0.0 


4.36 ± 


1.96 


SB 


9.49 ± 


4.25 


n 1^^+0-08 


o.i2l°o:S^ 


D4" 


3 


5.19 ± 1.73 


0.0 ± 


0.0 


400.0 ± 


0.0 


6.88 ± 


3.09 


SB 


14.97 ± 


6.71 


n 91+0-10 


o.2oi°:?i 



Table 3 

Properties of the CO emission of the EGNoG galaxies (described in Section 4.1). See Section B.2 for a description of the errors. 

" CO quantities are upper hmits. 



appropriate aco value, but we do not reflect this uncer- 
tainty in the error in the molecular gas mass. The error 
in Mmgas is calculated by propagating the error in L^q. 
The classification of galaxies using the specific star for- 
mation rate is described in more detail in Section 5.2. 
In Section 5.5, we discuss the use of a bimodal prescrip- 
tion for aco in this work and the the possibility of a 
continuously varying conversion factor. 

We report the molecular gas fraction and the molecular 
gas depletion time in the last two columns of Table 3. 
In keeping with previous work, we define the molecular 
gas fraction (/mgas) as the fraction of baryonic matter 
(ignoring atomic and ionized gas): 



/, 



M„ 



mgas 



M^+M, 



(5) 



mgas 



The molecular gas depletion time (rmgas) represents the 
time it would take to consume all the available molecular 



gas at the current star formation rate, defined: 



M 



'^dep 



mgas 



SFR 



(6) 



The errors in both of these quantities include both the 
error in the CO luminosity and the error in the stellar 
mass or SFR. 

The integrated CO emission maps of all detected galax- 
ies are presented in Figures 3-8. For each galaxy, the op- 
tical image (from SDSS) is in the left-most panel, the mo- 
ment (total intensity) map is in the left middle panel, 
the moment 1 (intensity-weighted mean velocity) map is 
in the right middle panel, and the integrated spectrum is 
in the far right panel. The dotted white ellipse shows the 
region in which pixels were summed (if not masked; see 
Section B.2 for details). For clarity, we have additionally 
masked the moment 1 map outside of the source region 
(dotted white ellipse) in Figures 3-8. The white bar 
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Figure 3. Integrated CO emission maps for detected bin A galaxies, part 1. The left panel shows the optical image with the galaxy name 
in the upper left. The moment and moment 1 maps are displayed in the left middle and right middle panels respectively The dotted 
white ellipse indicates the source region. In the moment map, the synthesized beam is indicated by the solid white ellipse in the lower 
left corner and a 10 kpc scale bar is given in the top right. The moment 1 map is masked outside of the source region ellipse for clarity. 
The far right panel shows the integrated spectrum of the galaxy: the solid blue line is calculated with masking, the dotted green line is 
without. 
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Figure 4. Integrated CO emission maps for detected bin A galaxies, part 2. Same as Figure 3. 
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Figure 5. Integrated CO emission maps for bin B galaxies, part 1. Same as Figure 3. 
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in the moment map shows the angular size of 10 kpc 
at the redshift of the galaxy. In the far right panel, the 
solid blue line shows the masked spectrum and the dot- 
ted green line shows the unmasked spectrum. Channels 
containing source emission are shaded yellow. 

The moment maps are produced using the 2(j clipped 
smooth mask method (see Section B.2). The moment 
map is a simple sum of the masked cube in the 'source' 
velocity channels. The moment 1 map is produced by 



summing the masked cube multiplied by the velocity in 
each channel, then normalizing by the moment value. 
Thus, the color in the moment 1 map indicates the 
intensity-weighted mean velocity in each pixel. The mo- 
ment 1 map is only calculated in pixels which have a 
positive moment value. Moment 1 pixels with values 
outside the velocity range covered by the source channels 
are masked. 
In most of the bin A galaxies, we spatially resolve the 
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Figure 6. Integrated CO emission maps for bin B galaxies, part 2. Same as Figure 3. 
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CO disk, illustrated by a velocity gradient across the 
emission region in the moment 1 map (indicative of a ro- 
tating disk). The morphology of the bin A CO emission 
is discussed further in Section 4.2. The bin B galaxies 
which were observed only in D configuration are spatially 
unresolved. However, galaxies Bl, B2, B3 and B7 were 
observed in C configuration as well, increasing the reso- 
lution by roughly a factor of 2. Of these four, galaxies 
Bl and B3 show evidence for ordered rotation but are 



only marginally resolved. Bin C galaxies are spatially 
unresolved in the CO (J = 1^0) line and marginally re- 
solved in the C0( J = 3^2) line. The C0( J = 3^2) 
moment 1 maps of the galaxies C2, C3 and C4 are sug- 
gestive of a rotating disk of gas. Bin D galaxies are not 
detected, but would be unresolved based on the sizes of 
the stellar disks. 

4.2. Morphology of the Low-Redshift Sample 
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Figure 7. Integrated CO emission maps for bin C galaxies in the CO (J = 1—^0) line. Same as Figure 3. 



For the EGNoG galaxies in redshift bin A, the spatial 
resolution of the CO observations allows us to comment 
on the morphology of the molecular gas disk at z = O.OS- 
O.l. The integrated CO intensity maps in Figures 3 and 
4 are suggestive of irregular morphology at first glance. 
However, due to the modest signal to noise ratio of these 
observations, it is difficult to infer the true shape of the 
underlying CO emission. 

To determine the effect of a modest signal to noise 
ratio on our CO emission maps, we produced a simple 
model disk galaxy at z = 0.08, inserted it into noise- 
only uv data channels from a bin A galaxy dataset and 
processed the data as we did the real data. We used 
a Milky- Way- like rotation curve, linear in the central 1 
kpc, flat outside. The distribution of the CO emission 
is exponential, with a scale radius of 3 kpc. We varied 
the total flux, rotation velocity of the flat part of the 
rotation curve (vflat), inclination angle {i) and physical 
radius of the molecular gas disk {Rqo). Models Ml to 
M4, shown in Figure 9, are representative of the results 



of varying all four parameters within reasonable ranges. 
The parameters used in the four models are listed in 
Table 4. 

The 'observed' model maps (right two panels in Figure 
9) suggest that at the signal to noise level typical of our 
bin A data, irregularities in the CO emission maps do not 
necessarily indicate an underlying irregular structure in 
the molecular gas disk. In most cases, our observations 
are consistent with the CO in the bin A galaxies being 
ordered in a rotating disk. We conclude that deeper ob- 
servations would be required to investigate the detailed 
morphology of the molecular gas disks in these systems. 

Two exceptions are galaxies A3 and A8, which exhibit 
disturbed morphologies that are not likely the result of 
low signal to noise ratios. Both galaxies show rotating 
molecular gas disks misaligned with the major axis of 
the optical disk, as well as CO emission outside the op- 
tical disk. These two galaxies are discussed in detail in 
Appendix E. 
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Figure 8. Integrated CO emission maps for (detected) bin C galaxies in the CO(J = 3 — )► 2) line. Same as Figure 3. 
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Table 4 

Model parameters: inclination angle (i), rotation velocity of the 

flat part of the velocity curve (-Uflat), radius of the CO disk (Rco) 

and total CO flux (S'co)- 



4.3. Non- Detections of the Bin D Galaxies 

We do not detect CO(J = 3 -^ 2) in any of the four 
bin D galaxies at z ~ 0.5. The Sa upper hmits on Sco 
are between 2.7 and 5.2 Jy km s~^, assuming a total ve- 
locity width of 400 km s~^. The spectra of the four bin 
D galaxies are shown in the top four panels of Figure 10 
at a velocity resolution of ^ 32 km s~^. The bottom 
panel shows the stacked spectrum, where each individ- 
ual spectrum is weighted by the inverse of the variance. 
Note that the spectra are stacked centered at the spec- 
troscopic redshifts given in Table 3, which have a typical 
uncertainty of 110 km s~^. The stacked spectrum does 
not indicate a detection; we calculate a 3(7 upper limit of 
1.5 Jy km s~^ for a 400 km s~^ wide signal. 

To place these hmits in context, we calculate the ex- 
pected CO flux of a galaxy with SFR ^ 70 M© yr"\ 
typical for the bin D galaxies. First we consider these 
galaxies to be normal, using a Milky way-like aco and 
rsi = 0.5. For Tmgas = 0.77 =b 0.57 Gyr (average from 
our literature sample; see Section 5.3), we expect a CO 



flux ^co = 1.2 — 8.4 Jy km s~^. Performing the same 
calculation using a starburst aco and the starburst av- 
erage Tmgas = 0.07 ± 0.04 Gyr (Section 5.3), we expect 
6'co = 0.76 - 2.7 Jy km s-^ We have used rsi = 0.5 
in this calculation but note that a larger value may be 
appropriate in a starburst system. We can further con- 
strain the expected CO flux using the typical stellar mass 
(log(M*) = 10.7) and the expected molecular gas frac- 
tion (see Section 5.4). For /mgas = 0.03 — 0.4, we ex- 
pect Sco = 0.13 — 2.6 Jy km s~^ for normal galaxies or 
*5'co = 0.5 — 11.0 Jy km s~^ for starburst galaxies. Note 
also that these estimates have an uncertainty of a factor 
of 1.5-2.5 from the uncertainty in the SFR and M* of the 
bin D galaxies, in addition to further uncertainty in the 
rsi value. 

Considering the range of expected CO fluxes, four non- 
detections at the depth of these observations is not un- 
reasonable. Given the large, overlapping ranges of ex- 
pected CO fluxes for starburst and normal galaxies de- 
tailed above, we do not draw any conclusions about the 
gas consumption properties of the bin D galaxies based 
on these upper limits. 

5. DISCUSSION 

We now look at the EGNoG galaxies in the context of 
starburst and normal galaxies at low and high redshift 
from the literature. We present the literature sample in 
the following section (Section 5.1), describe our identifi- 
cation and treatment of starburst and normal galaxies in 
the literature dataset (Section 5.2), then investigate the 
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Figure 9. Model galaxy maps. For each model (Ml - M4, parameters given in Table 4), the far left and left middle panels show the model 
galaxy moment and moment 1 maps respectively. In the right middle and far right panels are the moment and moment 1 maps from 
normal processing of the model galaxy inserted into noise-only data. 



molecular gas depletion times in normal and starburst 
galaxies (Section 5.3) and the evolution of the molecular 
gas fraction with redshift (Section 5.4). 

5.1. Literature Data 

The local galaxies of our literature data compilation in- 
clude normal spiral galaxies, LIRGs and ULIRGs. We in- 
clude 22 galaxies from Leroy et al. (2008), 10 of which are 
defined by the authors to be dwarf galaxies (M^ < 10^^ 
M0). Leroy et al. (2008) present CO data from the 
HERACLES and BIMA SONG surveys, stehar masses 
estimated from the SINGS survey if-band luminosities, 
and star formation rates calculated from a combina- 
tion of GALEX FUV and Spitzer 24 /im. We aug- 
ment this dataset with 47 galaxies (26 of which have 
M* < 10^^ M0) from the overlap of the Ha survey by 
Kennicutt et al. (2008) and the CO sample (taken from 
the literature) of Obreschkow & Rawhngs (2009). We 
estimate stellar masses, star formation rates and CO 
luminosities for these galaxies following Bothwell et al. 
(2009). For each galaxy, the stellar mass is calculated 
from the 5-band luminosity using a mass-to-light ratio 
estimated using the B-V color (Bell & de Jong 2001). 



The SFR is calculated from the Ha flux via the rela- 
tion given by Kennicutt (1998a), using an extinction cor- 
rection calculated from the B-hand luminosity following 
Lee et al. (2009). For more details, see Bothwell et al. 
(2009). 

We include a sample of 56 local LIRGs and ULIRGs, 
with CO data compiled from Gao & Solomon (2004), 
Solomon et al. (1997) and Sanders et al. (1991). We use 
stellar masses and SFRs from the following: the GOALS 
survey (Howell et al. 2010), with stellar masses calcu- 
lated from iT-band luminosties and SFRs estimated from 
the combination of IRAS infrared and GALEX FUV lu- 
minosities; da Cunha et al. (2010), with stellar masses 
and SFRs from SED fitting of mid-IR Spitzer/IRS spec- 
tra and multi-band photometry available from the liter- 
ature; and Hou et al. (2011) (table of values via private 
communication), with stellar masses from SED fitting 
of SDSS spectra. We estimate SFRs for the Hou et al. 
(2011) galaxies (for which only stellar masses are given) 
from the total infrared luminosity (Lir; 8-1000 /im) using 
SFR = 1.7 X 10-^^ Lfir (Kennicutt 1998a), where LpiR 
is the far- infrared luminosity in Lq) and i^m/i^FiR = 1.3 
(Gracia-Carpio et al. 2008). 
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Figure 10. Spectra of the four bin D sources within a circular 
region 2'^ in radius (top four panels), plus the stacked (average) 
spectrum, weighted by the inverse of the variance, in the bottom 
panel. 

At intermediate redshifts {z = 0.05 - 1), the EGNoG 
sample is augmented by the CO (J = 1^0) observations 
of 7 LIRGs (2 are upper hmits) at 2: ~ 0.4 by Geach et al. 
(2011), as weh as CO(J = 1^0), CO(J = 2^1) 
, and CO (J = 3 -^ 2) observations of 56 ULIRGs at 
0.2 < z < 1.0 by Combes et al. (2011) and Combes et al. 
(2013). We follow Combes et al. and use r2i = rs2 = 1, 
motivated by the fact that these are extremely star- 
forming systems and thus the gas is likely highly ex- 
cited. We exclude galaxies that are only observed in the 
CO(J = 4^3) hue due to the uncertainty of r^i. For 
the Geach et al. galaxies, stellar masses are from the K- 
band luminosities using a mass-to-light ratio determined 
from SED fitting using B VRIJK imdigmg. SFRs are from 
the total infrared luminosities, estimated from the PAH 
7.7 jam line (see Geach et al. 2009). Combes et al. (2011) 
estimate stellar masses for the 0.2 < z < 0.6 sample 
from K-corrected optical and near-IR luminosities. For 
the 0.6 < 2; < 1.0 sample. Combes et al. (2013) estimate 
stellar masses from SED-fitting of optical and near-IR lu- 
minosities from public catalogs. We estimate SFRs from 
the Combes et al. (2011, 2013) far-IR luminosities using 
the conversion of Kennicutt (1998b), as described above 
for a few of the LIRGs and ULIRGs. 

At high redshifts {z > 1), we include 67 SFGs pre- 
sented by Tacconi et al. (2012) in the IRAM Plateau 
de Bure High-z Blue Sequence CO (J = 3^2) Survey 
(PHIBSS) and the 16 of 26 SMGs from Bothweh et al. 
(2012) with stellar masses, at 2; < 3 (tables for both 
datasets via private communication). The SFGs were ob- 
served in the CO(J = 3^2) (Tacconi et al. 2010) and 
CO(J = 2^1) (Daddi et al. 2010; Magnelh et al. 2012) 



lines and their luminosities were converted to CO (J = 
1^0) luminosities using rsi =0.5 (Bauermeister et al. 
2013) and rsi = 0.8 (Dannerbauer et al. 2009). The 
SMGs from Bothwell et al. (2012) were observed in var- 
ious rotational transitions of the CO molecule ( Jupper = 
2-7) and their luminosities were converted to CO (J = 
1^0) luminosity using an excitation model calculated 
by Bothwell et al. from their entire sample, normal- 
ized to a common far-IR luminosity. The stellar masses 
for the Tacconi et al. (2012) SFGs are estimated from 
SED modeling and the SFRs are derived from a com- 
bination of the UV and IR luminosities or extinction- 
corrected Ha luminosity. Bothwell et al. (2012) use stel- 
lar masses from Hainline et al. (2011) (not available for 
all 26 galaxies), which come from stellar population syn- 
thesis fitting to observed-frame optical to mid-IR SEDs. 
We estimate SFRs for the SMC sample galaxies from 
the Bothwell et al. (2012) far-IR luminosities (which are 
calculated from the 1.4 GHz continuum utilizing the far- 
IR-radio correlation) using the conversion of Kennicutt 
(1998b), as described above for a few of the LIRGs and 
ULIRGs. 

While we have attempted to make this compilation of 
literature data as uniform as possible, the derived prop- 
erties are not all calculated in a consistent fashion. In 
converting CO luminosities to molecular gas masses, we 
have standardized the selection of aco value using the 
starburst classification of Equations 1 and 2 (described 
in the next section). The Jupper > 1 CO luminosities 
are converted to CO (J = 1 — > 0) luminosities using the 
ratios prescribed by the individual authors in the case of 
ULIRGs and SMGs. We use standard ratios for all other 
galaxies: rsi = 0.5 and r2i = 0.8. The M* and SFR 
values, however, are calculated by a variety of methods, 
as described above. This may have systematic effects on 
the analysis in the following sections. In the interest of 
comparing with as large a sample as possible, we include 
all the data described here in our analysis, but note that 
systematic effects may be present. 

5.2. Identifying Starburst Galaxies 

In Figure 11, we present the stellar mass and star for- 
mation properties of the literature sample and the EG- 
NoG survey galaxies. For the literature sample, the 
colors indicate approximate redshifts (2; ~ in blue, 
z ^ 0.05 — 0.5 in green, and 2; > 1 in red) and the 
symbols indicate galaxy type: normal galaxies (local spi- 
rals, Geach et al. 2011 galaxies at 2; = 0.4 and high- 
redshift SFGs) are marked with circles, local LIRGs with 
triangles, and starbursts (local ULIRGs, Combes et al. 
intermediate-redshift ULIRGs, and high-redshift SMGs) 
are marked with squares. The EGNoG galaxies are plot- 
ted with a unique symbol and color (magenta diamond) 
for emphasis. This color and symbol scheme will be used 
for the rest of this work. 

The left panel of Figure 11 shows the SFR versus M*, 
with the main sequence at three representative redshifts 
indicated by the blue (z = 0), green {z = 0.3), and red 
{z = 1.5) lines. Normal galaxies lie near the main se- 
quence curve at the corresponding redshift. In order 
to more precisely discern starburst galaxies from nor- 
mal galaxies, in the right panel of Figure 11 we plot the 
sSFR of each galaxy normalized by the main sequence 
sSFR (sSFRms; Equation 1) at the redshift and stellar 
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Figure 11. Star formation rate (left panel) and normalized specific star formation rate (right panel) versus stellar mass for our literature 
sample and the EGNoG galaxies. Colors indicate approximate redshifts and plotting symbols show the galaxy type (see text for full 
description). In the left panel, the main sequence at three representative redshifts is indicated by the blue, green and red lines. In the 
right panel, the sSFR is normalized by the main sequence sSFR (for the relevant z and M*). Starburst (SB) galaxies lie above the higher 
dashed line. 



mass of the galaxy. The horizontal black line indicates 
the main sequence relation, with the dashed lines show- 
ing a spread of 0.6 dex. Starburst galaxies lie above the 
higher dashed line, with sSFRs larger than 4 times the 
main sequence sSFR (as defined in Equation 2). In most 
cases, our starburst classification agrees well with expec- 
tations (e.g. ULIRGs are mostly classified as starbursts), 
with the exception of the SMG population, which we 
find to be mostly consistent with the main sequence de- 
fined in Equations 1 and 2. However, we note that since 
most of the SMGs included here are at z > 2, this re- 
sult depends sensitively on the high-redshift behavior of 
the sSFRms(^) prescription, which is more uncertain and 
may fiatten at these redshifts (see discussion in Section 
2.1). 

Using this starburst criterion, we standardize the con- 
version of CO luminosity to molecular gas mass in our 
literature sample, as we did for the EGNoG galaxies in 
Section 4.1. For galaxies that are classified as starburst 
based on their specific star formation rates, we use the 
value observed in starburst systems: Qfco(ULIRG) = 0.8 
Mo (K km s-i pc2)-i (local ULIRGs: Scoville et al. 
1997, Downes & Solomon 1998; high-redshift SMGs: 
Tacconi et al. 2008). For all other galaxies, we use a 
Milky Way-like conversion factor, aco (Milky Way) = 3.2 
M© (K km s-^ pc^)"^ (e.g. Dame et al. 2001). These 
conversion factors do not explicitly include He; see Equa- 
tion 4. This bimodal prescription for the conversion fac- 
tor is discussed further in Section 5.5. 

5.3. Gas Depletion Time 

Using our compilation of data, we investigate the 
molecular gas depletion time (rmgas) for starburst and 



normal galaxies. Excluding dwarf galaxies, the aver- 
age molecular gas depletion time for the normal (non- 
starburst) galaxies is 0.76±0.54 Gyr. The average for the 
starburst galaxies is 0.06±0.04 Gyr. (In calculating these 
averages, we have excluded galaxy Q2343-MD59 from 
Tacconi et al. 2012, as it is an obvious outlier.) These 
averages generally agree with previous work (see refer- 
ences for the literature sample in Section 5.1). However, 
we note that a depletion time calculated from molecu- 
lar gas and SFR surface densities may differ from ours 
(calculated from the total molecular gas mass and SFR) 
due to different scale lengths of the CO and SFR dis- 
tributions and trends in gas depletion time with galaxy 
size (larger galaxies will dominate an area-weighted av- 
erage gas depletion time). This is the case for the sample 
of Leroy et al. (2008), in which the authors calculate an 
average molecular gas depletion time of 1.9 Gyr using 
surface densities, while the average depletion time using 
total quantities would be 1.3 Gyr (which lies within the 
error of our average value for normal galaxies). 

In the left panel of Figure 12, we plot CO luminos- 
ity versus SFR to illustrate that starburst and normal 
galaxies fill the same space in this plane, making the dif- 
ferentiation of starburst from normal galaxies impossible 
using only CO luminosity and SFR. In addition to the lit- 
erature sample described in Section 5.1, we include here 
local spirals, LIRGs and ULIRGs from Yao et al. (2003), 
Mao et al. (2010) and Papadopoulos et al. (2012) as gray 
circles, triangles and squares respectively. These datasets 
do not include stellar masses and thus are only included 
in this plot for illustrative purposes (we estimate SFRs 
from the far-IR luminosities using the Kennicutt 1998b 
conversion). The dashed lines indicate the area of the 
plane occupied by normal galaxies using aco (Milky Way) 



The EGNoG Survey 



17 



11 



10 



E 







^.oi■■■<^ 



111 ^ 



typical SFR error 





local spirals 




local LIRGs 




local ULIRGs 




EGNoG 




Geach2011 




Combes2013 




high-z SFGs 




high-z SMGs 



10.5 



10.0 



normal 
starburst 



typical SFR error 



-10 12 

log(SFR/MQ yrM 




12 

log(SFR/MQ yrM 



Figure 12. CO luminosity (^qq' ^^^^ panel) and molecular gas depletion time (rmgas; right panel) versus SFR. In the left panel, galaxies 
from Yao et al. (2003), Mao et al. (2010), and Papadopoulos et al. (2012) are plotted (gray points) in addition to the literature sample. 
Dashed (dotted) lines show the expected region in L^q-SFR space occupied by normal (starburst) galaxies. In the right panel, the horizontal 
solid (dot-dashed) line indicates the average Tmgas value for normal (starburst) galaxies, the gray region indicating ztlcr. Starburst galaxies 
are plotted as un-filled symbols in the right panel. 



and the average Tmgas for normal galaxies given above. 
The dotted lines indicate the area occupied by starburst 
galaxies (using aco(ULIRG) and the starburst Tmgas), 
which overlaps the normal galaxy area significantly. This 
overlap occurs because the differences in aco and Tmgas 
between starburst and normal galaxies cancel each other 
out in this plane: as the effective aco increases, a galaxy 
moves toward the lower right and as the effective Tmgas 
increases, a galaxy moves toward the upper left. 

In the right panel of Figure 12, the molecular gas de- 
pletion time is plotted against SFR, with the average 
value for normal (starburst) galaxies indicated by the 
soHd (dot-dashed) line, la error bars on the averages 
are indicated by the gray shaded regions. The CO to 
H2 conversion factor used for each data point is indi- 
cated by the fill of the symbol: filled symbols are normal 
(ofco (Milky Way)) and un-filled symbols are starburst 
(Qfco(ULIRG)) according to the sSFR criterion in Equa- 
tion 2 (plotted in the right panel of 11. Only galaxies 
that can be classified this way (galaxies for which the 
M* and SFR have been estimated) are included in this 
plot. Note that the Tmgas error bars plotted only repre- 
sent the error associated with the CO measurement and 
not the SFR error since we do not have errors on all the 
SFR measurements in the literature sample. We estimate 
that the typical SFR error is ±50%. 

The population of non-starbursting, local galaxies in 
the lower left corner of the right panel are dwarfs (M* < 
10^^), for which a larger conversion factor may be more 
appropriate due to decreasing metallicity. Leroy et al. 
(2011) find that the conversion factor begins to in- 
crease significantly as the metallicity falls below 12 + 
log(0/H) = 8.2-8.4, which is expected for M* < 10^ M© 



(Tremonti et al. 2004). This change in conversion fac- 
tor would shift the dwarf galaxy points to higher Tmgas • 
Therefore, we do not include dwarf galaxies in the calcu- 
lation of the average Tmgas values. 

5.4. Evolution of the Molecular Gas Fraction 

In order to calculate the expected evolution of the 
molecular gas fraction, we first consider the ratio of the 
molecular gas mass to the stellar mass (rmgas), which can 
be written as the product of the molecular gas depletion 
time and the specific star formation rate: 



M 



mgas 



' mgas 



M* 



^mgas SFR 

SFR M* 



Tmgas X sSFR (7) 



In Section 2.1, we motivated and defined the sSFR of 
the main sequence as a function of redshift and stellar 
mass. As illustrated in the right panel of Figure 11, main 
sequence galaxies lie roughly within a factor of 4 (±0.6 
dex) of the main sequence relation and starburst galaxies 
lie at sSFRs 4 to 30 times larger than the main sequence 
value. In the previous section, we found an average Tmgas 
of 0.76 ±0.54 Gyr for normal galaxies and 0.06 ±0.04 Gyr 
for starburst galaxies. Using these typical Tmgas values 
and sSFR ranges (relative to the sSFRms; see Equation 
1), we get: 



^mgas (normal) 



(0.05 ± 0.04)7/norm(l ± zf-^M'^^^ 



^mgas (starburst) = 

(0.04 ± 0.03) (^)(1±^)^-"M^ 



-0.2 

*,ii 
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Figure 13. Molecular gas ratio (rmgas) and molecular gas fraction (/mgas) versus redshift. We plot log{z + 0.1) on the x-axis. The 
average value from the COLDGASS survey is plotted as a cyan upside-down triangle. Starburst galaxies are plotted as un-filled symbols. 
Average values (with la errors) for each of 6 redshift bins are plotted as un-filled black upside-down triangles. The shaded gray area shows 
the expected behavior for normal galaxies, with the solid black curve indicating the average. The dotted black lines bound the expected 
behavior for starburst galaxies. 



where M^^u = M^/{10^^ M©) and r^norm and t^sb are the 
ranges of sSFR/sSFRms values appropriate for normal 
and starburst galaxies, respectively: r^norm = 1/4 to 4 
and T^sB = 4 to 30. In each of Equations 8 and 9, the 
scatter in Tmgas is reflected in the error of the flrst term 
(zbO.04 and ±0.03). The corresponding molecular gas 
fraction may be calculated from rmgas- /mgas = (1 + 



' mgas 
'"mgas 



) . The result is that the expected behavior of 
and /mgas is remarkably similar for normal and 
starburst galaxies. 

Figure 13 shows the data and expectation for rmgas 
(left panel) and /mgas (right panel) versus redshift, plot- 
ting log{z + 0.1) on the x-axis in order to spread out the 
low to intermediate redshifts for clarity. As in the right 
panel of Figure 12, the starburst galaxies (identifled us- 
ing Equation 2) are plotted as un-filled symbols. Note 
that the error bars plotted only represent the error as- 
sociated with the CO measurement (and not the stellar 
mass error), since we do not have errors on all the stellar 
mass measurements in the literature sample. We esti- 
mate that a typical stellar mass error is 0.1-0.2 dex. We 
also plot the average value at z = 0.025 — 0.05 reported 
by the COLDGASS survey (Saintonge et al. 2011) as a 
cyan upside-down triangle. We plot the expected rmgas 
and /mgas versus redshift for normal galaxies (Equation 
8) as the solid black curve, with the gray shaded area 
indicating the range of possible values (we take a nom- 
inal M* of 10^^ M0). The range for starburst galaxies 
(Equation 9), indicated by the dotted black curves, over- 
laps the expected range for normal galaxies. 

The data points are binned into 7 redshift ranges (num- 
ber of detected galaxies in parenthesis): z = — 0.05 
(110), 0.05 - 0.1 (19), 0.15 - 0.25 (18), 0.25 - 0.45 (18), 
0.45 - 1.0 (15), 1.0 - 1.7 (53), and 1.9 - 2.7 (28). The 
average values (for all detections) are plotted as un-filled 



black upside-down triangles, with vertical lines showing 
the standard deviations (for the z = 0.25 — 0.45 bin, the 
obvious outlier from Combes et al. 2011 at rmgas = 3.8 
is excluded). All redshift bin averages lie within the ex- 
pected range for both normal (gray shaded region) and 
starburst (dotted lines) galaxies, and the average expec- 
tation curve (rmgas = sSFRms x (0.76 Gyr)) lies within 
the dispersion (vertical black lines). 

We note that the three intermediate redshift bins {z = 
0.05-0.01, 0.15-0.25, and 0.25-0.45) lie systematically 
higher than the average expectation curve. However, the 
number of galaxies in each of these bins is small and is 
dominated by SFGs (EGNoG and Geach et al. galaxies), 
which we have shown to lie close to the starburst cutoff 
(see Figure 11). Therefore, these points are sensitive 
to the exact starburst criterion used (which determines 
aco)' For example, if we consider galaxies with sSFR 
> 2 sSFRms to be starbursts, the three intermediate 
redshift bin average values decrease by ~ 0.2 dex, falling 
more in line with the average expectation curve. 

In summary, while the EGNoG galaxy redshift bins lie 
systematically higher than the average expectation curve 
for the starburst criterion (and resulting aco choices) 
adopted in this work, a reasonable variation of this cri- 
terion results in a decrease of ^ 0.2 dex, bringing the 
averages in line with expectation. We also note, but do 
not investigate, that systematic effects may be present 
due to the use of different methods for calculating M* 
and SFR, and from instrumental sensitivity limiting CO 
detections to gas-rich galaxies at intermediate and high 
redshift. Overall, we conclude that the data agree well 
with the behavior of rmgas and /mgas expected from a 
simple prescription of sSFR and Tmgas with z and M* in 
star- forming galaxies. While the EGNoG bin D galaxies 
(non-detections) have not been included in the redshift 
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bin averages discussed above, we note that the four upper 
hmits are included in Figure 13 and agree with expecta- 
tions as weh. 

5.5. A Bimodal Conversion Factor Prescription 

In this work, we have used a bimodal prescription for 
the conversion factor aco in normal and starburst galax- 
ies. We have included some local low-mass galaxies in 
our literature sample, but we do not include these sys- 
tems in the calculation of average values since a differ- 
ent aco is likely appropriate due to decreasing metal- 
licity (Leroy et al. 2011; Narayanan et al. 2012). While 
our bimodal prescription (excluding low-mass galaxies) 
likely does not describe the true conversion factor in all 
of these galaxies (e.g. simulations by Narayanan et al. 
2012 suggest a continuous variation of aco from normal 
to starburst galaxies), it appears to capture the typical 
behavior. A recent study by Magnelli et al. (2012) using 
observed dust masses to calculate the conversion factor 
in z > 1 star-forming galaxies found smaller conversion 
factors in starburst galaxies (classified by sSFR, as in 
this work), but it was unable to distinguish between a 
step function and a continuous variation of aco as a 
function of sSFR/sSFRms- In this work, we have used 
only total measurements of the galaxy properties, and 
we note that with resolved CO measurements, it would 
be possible to apply a value of aco more tailored to lo- 
cal conditions (like the prescriptions of Narayanan et al. 
2012; Bolatto et al. 2013), which could tighten the ob- 
served relations. 

Independent of this issue of bimodal or continuous 
conversion factor is the question of whether the lo- 
cal values for normal and starburst conversion factors 
may be extended to high redshift. Both simulations by 
Narayanan et al. (2012) and observations by Daddi et al. 
(2010) suggest that the conversion factor is lower than 
the Milky Way value in normal star-forming galaxies at 
high redshift. We cannot answer this question with the 
current dataset, but note that using a reduced aco in 
high-redshift SFGs or a continuous variation of aco from 
normal to starburst galaxies would not change the gen- 
eral conclusions. 

6. CONCLUSIONS 

In this paper, we present the EGNoG survey: CO ob- 
servations of 31 galaxies from 2: = 0.05 to 2: = 0.5. We de- 
tect 24 of the 31 observed galaxies, providing integrated 
CO maps, spectra, and luminosities for detections and 
3(7 upper limits for non-detections. 

We place the EGNoG galaxies in context with a sam- 
ple of 185 normal and starburst galaxies at low and high 
redshift collected from the literature. We standardize the 
comparison of the EGNoG and literature galaxies using 
a simple prescription for molecular gas mass calculation. 
Each galaxy is classified as normal or starburst based 
on its sSFR relative to the sSFR of the main sequence 
of star-forming galaxies at that redshift and stellar mass 
(see Equations 1 and 2). We then use a bimodal pre- 
scription for aco based on this classification (a Milky 
Way- like value for normal galaxies, a ULIRG-like value 
for starbursts), and calculate the molecular gas mass, 
molecular gas depletion time, and molecular gas fraction 
of each galaxy. 



Using this prescription, we find an average molecular 
gas depletion time of 0.76 ± 0.54 Gyr for normal galax- 
ies and 0.06 ± 0.04 Gyr for starburst galaxies. We cal- 
culate an average molecular gas fraction of 10-20% at 
the intermediate redshifts probed by the EGNoG sur- 
vey {z = 0.05 — 0.5). By expressing the molecular gas 
fraction in terms of the sSFR and molecular gas deple- 
tion time (/mgas = (1 + [^mgas X sSFR]"^)"^), and using 
typical ranges of sSFR and Tmgas for starburst and nor- 
mal galaxies, we calculate the expected evolution of the 
molecular gas fraction with redshift. We find that the 
expected behavior for normal and starburst galaxies is 
remarkably similar and agrees well with the EGNoG and 
literature data from z = out to z ^ 2.5. 
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APPENDIX 

A. STELLAR MASS AND SFR PROBABILITY DISTRIBUTION FUNCTIONS 

For each of the EGNoG galaxies in redshift bins A-C, we show the stehar mass and SFR probabihty distribution 
functions in Figure Al in the top and bottom panels respectively. These distributions come from the MPA-JHU group 
(see text) for SDSS DR7 galaxies. Data points give the 2.5, 16, 50, 84 and 97.5 percentiles of the distribution. The 
median, mean and mode are indicated by the vertical solid, dashed and dotted lines respectively. For galaxies Al, 
B4 and CI, a duplicate SDSS source exists (as a result of SDSS automated source-finding). The duplicate source is 
plotted in red in Figure Al. For each of the quantities plotted in Figure Al (stellar mass in the top panels, SFR in 
the bottom panels), the x-axis range is constant. 
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Figure Al. Stellar mass (top) and SFR (bottom) probability distribution functions of EGNoG galaxies selected from SDSS (bins A, B 
and C). For galaxies Al, B4 and CI, a duplicate SDSS source exists, over-plotted in red. 
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Calibrator 


Obs. 


Freq. 




Flux 


Primary Calibrator 




Name 


Band 


(GHz) 


Date Range 


(Jy) 


(Flux in Jy) 




UUU6-U63 


3mm 


105 


Oct 2010 


1.6 


MWC349 (1.28 Jy) 


0010+109 


3mm 


109 


Nov 2010 


0.2 


Uranus 


0108+015 


3mm 


97 


Sep 2011 


2.0 


Uranus; Neptune 




0108+015 


3mm 


97 


Feb 2012 


1.7 


Uranus 




0854+201P 


3mm 


107 


Nov 2010 


5.2 


3C84 (11 Jy) 




0854+201^ 


3mm 


98 


Aug- Oct 2011 


4.0 


Mars 




0854+201P 


3mm 


98 


Nov 2011 


5.2 


Mars 




0854+201P 


3mm 


98 


Feb 2012 


5.2 


Mars 




0854+201 


1mm 


258 


Aug 2011 


2.0 


Mars 




0854+201 


1mm 


254 


Apr 2012 


4.0 


Mars 




0920+446 


3mm 


106 


Nov 2010 


1.9 


3C84 (11 Jy) 




0920+446 


3mm 


106 


Feb 2012 


1.0 


3C84 (15 Jy) 




0958+655 


3mm 


105 


Nov 2010 


1.1 


3C273 (12 Jy) 




0958+655 


3mm 


105 


Feb 2012 


1.7 


MWC349 (1.28 Jy) 




1058+015^ 


3mm 


107 


May 2011 


4.6 


3C84 (10 Jy) 




1058+015^ 


3mm 


99 


Nov 2011 


4.1 


Mars 




1058+015^ 


3mm 


99 


Feb 2012 


3.6 


Mars 




1058+015^ 


3mm 


99 


Apr 2012 


2.8 


Mars 




1058+015 


1mm 


266 


Aug 2011 


2.6 


Mars 




1058+015 


1mm 


226 


Apr 2012 


2.0 


Mars 




1159+292 


3mm 


105 


Apr-May 2011 


1.4 


3C84 (10 Jy); 3C273 (10 Jy) 




1159+292 


3mm 


105 


Feb 2012 


0.7 


Mars 




1224+213 


1mm 


252 


Apr 2012 


0.6 


Mars 


1310+323 


3mm 


90 


Sep-Nov 2011 


1.7 


Mars 




1310+323 


3mm 


98 


Apr 2012 


1.4 


Mars 




1310+323 


1mm 


266 


Aug 2011 


0.6 


MWC349 (2.0 Jy) 




1357+193 


3mm 


105 


Apr-May 2011 


0.7 


MWC349 (1.28 Jy) 




1357+193 


3mm 


105 


Sep-Nov 2011 


0.8 


Mars; MWC349 (1.28 Jy) 




2134-018 


3mm 


106 


Apr 2011 


1.2 


Uranus 


2232+117 


3mm 


106 


Apr 2011 


1.3 


MWC349 (1.28 Jy) 




2232+117 


3mm 


97 


Aug- Sep 2011 


1.5 


Uranus; MWC349 (1.23 Jy) 




2232+117 


3mm 


97 


Feb 2012 


1.7 


Uranus; Neptune; MWC349 (1.23 


Jy) 


3C273 


3mm 


98 


Apr 2012 


6.8 


Mars 


3C454.3 


3mm 


107 


Oct 2010 


39.5 


Uranus 



Table Bl 

Summary of calibrator fluxes used to set the flux scale of EGNoG data. For each flux, the observing band (1 or 3 mm), frequency, and 
date range are listed. The last column lists the primary calibrator used to set the flux. For non-planet calibrators, the flux used (at the 

frequency given in the third column) is given in parentheses. 
P Indicates the calibrator is linearly polarized, requiring additional reduction steps (only applicable for 3 mm data) detailed in 

Appendix C. 

B. DATA REDUCTION AND FLUX MEASUREMENT 
B.l. Data Reduction 

Each dataset was reduced and calibrated as follows. The data were flagged for antenna - antenna shadowing and 
any other issues present during the observation (e.g. high system temperature, high gain, tracking problem, etc.). The 
instrument bandpass was calibrated with mf cal on a bright passband calibrator. The time-dependent antenna gains 
(from atmospheric variation) were derived by performing a self cal on the phase calibrator with an averaging interval 
of 18 minutes (the timescale of switching between the source and phase calibrator). For galaxies A5, All, B3, B4, B7, 
CI and C2, the phase calibrator (0854+201 or 1058+015) was up to 10% polarized, which required additional steps 
in the reduction of the 3mm data (observed with linearly polarized feeds). This is described in detail in Appendix C. 

For each dataset, the flux of the phase calibrator was set during the antenna gain calibration in order to properly 
set the flux scale of the data. The flux of each phase calibrator was assumed to be constant over timescales of weeks, 
and was therefore determined from the best datasets of the survey using bootf lux on bandpass-calibrated, phase-only 
gain-calibrated data using a planet (Uranus, Neptune or Mars) or MWC349 as a primary flux calibrator. In some 
cases, none of these flux calibrators was available so a bright quasar (3C84 or 3C273) was used instead, with the flux 
estimated using historical flux monitoring data at CARMA^. Table Bl provides the flux used for each phase calibrator 
and the primary calibrator used. For non-planet primary calibrators, the flux used is given in parenthesis. The flux 
of MWC349 was calculated assuming 1.2 Jy at 92 GHz with a spectral index of 0.5, the typical value from historical 
flux monitoring at CARMA^'^. The brightness temperature (T5) of the planets was set by the CARMA system. The 
brightness temperature of Mars comes from the Caltech thermal model of Mars (courtesy of Mark Gurwell), which 
includes seasonal variations in temperature. This model can be accessed in MIRIAD using marstb. For Uranus and 



^ See CARMA Memo #59 (Bauermeister et al. 2012) for a de- 
scription of flux monitoring at CARMA 



^ A history of MWC349 fluxes mea- 

sured at CARMA is maintained at 

http : //cedarf lat . mmarray . org/f luxcal/primary_sp_index . htm 
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Neptune, the CARMA system uses the fohowing power laws: 

T.(Uranus) = 134.7 (^^^)-°-"^K (Bl) 

/ V \ -0.350 

Tb(Neptune) = 129.8 [ ^^^ ^^J K (B2) 

Data cubes were produced combining ah fuhy cahbrated datasets of each source. We used invert, weighting the 
visibihties by the system temperature as weh as using a Briggs' visibihty weighting robustness parameter (Briggs 
1995) of 0.5. Since CARMA is an inhomogeneous array (these data use both the 10 and 6 m dishes), we also use 
options=inosaic in the invert step to properly handle the three different primary beams patterns (lOm-lOm, 10m- 
6m and 6m-6m). All observations consisted of a single pointing. The resulting cubes are primary-beam-corrected. If 
channel averaging was required (based on the strength of the CO line), this was done in the invert step. 

We deconvolved each image with mossdi (the mosaic version of clean), cleaning down to the rms noise within 
a single channel, within a cleaning box selected by eye to include only source emission. We cleaned only channels 
containing visible source emission. Cleaning down to a specified noise level is preferred to using a set number of clean 
iterations due to the nature of the spectral line emission: some channels will contain more flux and therefore require 
more clean iterations. In tests using a model source of known flux inserted into real data (emission- free channels), we 
found a Icr cutoff to best extract the true source emission without overestimating the flux over a range of detection 
significance levels, with a 10-30% uncertainty in the recovered flux (depending on the significance of the signal). The 
final cleaned cubes were produced by rest or, which convolves the clean component cube with a clean beam (calculated 
by fitting a Gaussian to the combined mosaic beam given by mospsf ) and adds the residuals from the cleaning process. 

B.2. Flux Estimation 

Total source fluxes in the CO hues are calculated by summing 'source' pixels in each 'source' velocity plane of the 
final, calibrated, cleaned cubes. The source velocity planes are selected by eye. The source pixels are those within 
the 'source' region which are not masked by our smooth mask, adapted from the masked moment method of Dame 
(2011). The smooth mask is created by applying a 2cr clip to a smoothed version of the image: Hanning smoothing 
is done along the velocity axis, and each velocity plane is convolved spatially with a Gaussian beam twice the size of 
the original synthesized beam. This smoothed mask is used to exclude noise pixels but still capture low-level emission 
that a simple a-clip (without smoothing) would miss. In our own testing (using a model source of known flux inserted 
into real data), we found a 2cr clip to best reproduce the true flux over a range of detection significance levels. 

An appropriate source region size was selected to recover all of the flux without including the negative bowl. Since 
we do not have single-dish data to complement the interferometric data presented here, our datasets are missing the 
shortest uv spacings. As a result, the emission in the clean maps sits in a negative bowl, which will affect the measured 
flux (calculated by summing pixels within a given radius). The bin A galaxies are observed with a sufficiently small 
beam so that the galaxies are resolved. In this case, we use an elliptical source region selected by eye to enclose 
the spatial extent of the CO emission of each galaxy. In the case of the bins B and C galaxies, the molecular gas is 
unresolved, so we use a circular source region with radius set to the radius at which the radial profile of the enclosed 
flux first peaks, thereby excluding the negative bowl. The radial profiles of the enclosed flux for bin B and C sources 
are shown in Figure Bl. The negative bowl is most evident in the unmasked profiles (dotted lines in Figure Bl), in 
which the enclosed flux peaks and then decreases with radius. We have selected the following source region radii by 
eye (approximately the radius at which the enclosed flux first peaks): 4.5'' for bin C sources in the C0( J = 3^2) line 
(1 mm data); b" for sources Bl, B2, B3, B7 and B9; Q" for B4 and B8; 6.5'' for bin C sources in the CO(J = 1^0) 
line (3 mm data); 7" for B6 and BIO; 8" for B5. The source region radius of each galaxy in bins B and C is indicated 
by the vertical dashed line in each panel of Figure Bl. For all galaxies, the source region is indicated by the dotted 
white ellipse in the moment and moment 1 maps of Figures 3-8. 

The error in each flux measurement is estimated from the standard deviation of the measured fluxes using different 
velocity channel averaging and starting channels, using three different methods of calculating the flux in each case. 
The three methods are: the 2cr masking technique described above, the same masking technique with a 3(7 clip, and 
the simple addition of all pixels (no mask) within the source region. 

We performed extensive testing of our analysis technique to choose the parameters of the reduction to eliminate 
systematic offsets and minimize the uncertainty due to noise (as described above). We find a 10-30% error in the flux 
measurement in each channel coming from noise in the data and the adopted reduction and analysis steps. From this, 
we adopt an average uncertainty of 20% in the flux estimated in each channel, which results in an uncertainty in the 
total flux of 2^%N~^'^ ^ where N^h is the number of velocity channels in which the flux is summed. For the total flux 
values reported in Table 3, we estimated the errors from the variations in the fluxes calculated using different channel 
averaging, flux measurement methods, etc. (as described above), consistent with the 20%(A/'ch)~^"^ expected. 

Further, these data suffer from systematic effects due to the absolute flux calibration and primary beam correction. 
We set the flux scale in our dataset using a primary flux calibrator, the flux of which is only known to ^ 20%. In 
the primary beam correction of the dataset, pointing and focus errors at the time of the observations as well as errors 
in the primary beam model can significantly reduce image fidelity, leading to errors in the measured fluxes of ^ 20% 
(see Square Kilometer Array (SKA) Memo #103, Wright & Corder 2008). Combining these errors in quadrature, we 
estimate that our flux measurements suffer from additional uncertainties of up to ^ 30%. We consider all these factors 
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Figure Bl. Enclosed flux versus radius for bin B and C sources. The masked (unmasked) data are shown by the sohd (dotted) hne. The 
vertical dashed line indicates the source region radius adopted for each galaxy. For galaxies C2, C3 and C4, the CO (J = 1^0) line data 
are shown in black and the CO (J = 3 — )► 2) line data in blue. 



in the presentation of our data in Table 3: for the hne flux (5'co), we report the measured error; for L'qq^ we include 
an additional 30% error (added in quadrature to the measured error in ^co)- 

C. POLARIZED CALIBRATORS 

In this section, we describe how we calibrated datasets which used a strongly polarized source for gain calibration. 
This applies to the 3mm data only, which were taken using single-polarization, linearly polarized receivers. The 1mm 
data used dual-polarization, circularly polarized feeds and are thus not affected by the strong linear polarization present 
in some calibrators. 

For Stokes I intensity /, linear polarization fraction p^^, and polarization angle PA, the observed YY intensity as a 
function of parallactic angle x is described by 



YY = I{l-pquCo^[2{PA-x)]} 



(Cl) 



Thus, the observed YY amplitude varies sinusoidally with parallactic angle in a dataset. Any amplitude self-calibration 
on the gain calibrator would assume constant amplitude with time and yield higher or lower gains according to this 
amplitude variation. Since this effect can be as large as 10% in some cases, we corrected for it in the following way. 
After a phase-only self cal, we used a modifled version of uveal to remove the known polarization effect from the 
calibrator data and then performed an amplitude self cal. The resulting gains should represent the atmospheric 
variations during the track, with the influence of the polarization of the calibrator removed. 

This procedure was used for two calibrators: 0854+201 and 1058+015. We do not assume stability over long 
timescales in the polarization parameters of the calibrators (and in fact flnd evidence for signiflcant variation with 
time), so the polarization parameters for each source are derived in each epoch from the best datasets. The flux scale 
in most of the datasets was calibrated using Mars, with brightness temperatures from the Caltech thermal model by 
Mark Gurwell, which includes seasonal variations. Datasets from November 2010 and May 2011 use 3C84 as the flux 
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Figure CI. Polarization of calibrators 0854+201 (left) and 1058+015 (right). In each panel, flux versus parallactic angle is plotted for the 
best datasets at each epoch (indicated in the upper left). The solid black line shows the best flt, with the dotted black line indicating the 
best-flt Stokes I flux. Calibrator 0854+201 (left) was observed at four epochs: at 107 GHz in November 2010 and at 98 GHz in September 
2011, November 2011 and February 2012. Calibrator 1058+015 (right) was observed at four epochs as well: at 107 GHz in May 2011 and 
at 99 GHz November 2011, February 2012 and April 2012. 

calibrator. Based on flux monitoring at CARMA'^, we assumed a flux of 11 Jy in November 2010 and 10 Jy in May 
2011 (both at 107 GHz). 

The fits in the different epochs are shown in Figure CI, with 0854+201 on the left and 1058+015 on the right. Flux 
versus parallactic angle is plotted, with data points for each dataset marked with a different symbol (indicated in 
the legend) and the best fit line (Equation CI) in solid black. The dotted black line shows the best fit Stokes I flux. 
The observing frequency and best-fit parameters are given in the upper right corner of each panel. The polarization 
parameters used for each calibrator, in each epoch, are given in Table CI. 

While the fit for 0854+201 in November 2010 is not well constrained due to limited parallactic angle coverage, it 
agrees with subsequent epochs, so we use it. The fit for 0854+201 in September 2011 is not very good due to noisy 
data, but the good fit from November 2011 (scaled to the average flux in September 2011; black dashed line) appears 
to match the September 2011 data well. Therefore, we use the November 2011 polarization fraction and angle for 
September 2011, with an average Stokes I flux of 4.0 Jy. 

The data for 1058+015 from November 2011 to April 2012 suggest a radical change in polarization, from 9% 
polarized at 19 degrees to 8% polarized at 94 degrees over the course of 6 months. The data from February 2012 show 
no indication of polarization (in fact, the best fit is a negative polarization fraction). This disappearance and (rotated) 



'^ See CARMA Memo #59 (Bauermeister et al. 2012) for a de- 



scription of flux monitoring at CARMA 



26 



Bauermeister et al. 



Calibrator 


Epoch 


^obs 


Stokes 1 


Pqu 


Angle 






(GHz) 


Flux (Jy) 




(deg) 


0854+201 


Nov 2010 


106.6 


5.2 


0.10 


70 




Sep 2011 


98.1 


4.0 


0.10 


67 




Nov 2011 


98.0 


5.2 


0.10 


67 




Feb 2012 


98.0 


5.2 


0.08 


54 


1058+015 


May 2011 


107.0 


4.6 


0.08 


19 




Nov 2011 


98.9 


4.1 


0.09 


19 




Feb 2012 


98.9 


3.6 


- 


- 




Apr 2012 


98.3 


2.8 


0.08 


94 



Table CI 

Polarized calibrator parameters. The values given in the table are those that were used in the reduction of the relevant datasets. In some 

cases, these are not the same as the best fit values. See text for discussion. 



reappearance of the linear polarization component of 1058+015 is coupled with a significant decrease in flux. 

D. f/y-SFECTRA 

We show the ix'^-spectra for all detected galaxies in Figures Dl and D2. For each galaxy (indicated in the top left), 
the top (bottom) panel shows the vector-averaged Real and Imaginary amplitudes (phase) of the calibrated uv data 
versus velocity. In Figure D2, the CO(J = 1^0) hue is shown in the top panels and the CO(J = 3^2) line in the 
bottom panels. For galaxy A8 and the CO(J = 3^2) line in galaxies CI and C2, the spectrum is calculated at a 
position offset from the SDSS position in order to coincide with the observed peak of the CO emission. For a compact 
source at the center of the field of view, the Real part shows the signal without a noise bias, and the Imaginary part 
shows the noise. In all cases (except CO(J = 3 -^ 2) in galaxy CI), we see coherent emission (larger Real amplitudes 
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Figure Dl. u-u-spectra of bin A and B galaxies. The top panel shows the vector- aver aged Real (solid) and Imaginary (dotted) amplitudes 
(Real,Imag in mJy beam"-*^) and the bottom panel shows phase (ph, degrees) of the calibrated uv data versus velocity (km s"-*^). For 
galaxy A8, the spectrum is calculated at a —0.8^^, 0.8^^ offset from the SDSS position (indicated in the upper right corner of the A8 panel). 



The EGNoG Survey 



27 



> 10 
F 8 



oi 30 
I 20 



-10 



S 

°- -100 



CI 



CO(l-O) 



CO(3-2) 




C2 



j^ 



CO(l-O) 



»■•■••■» • •■•■■ 



CO(3-2) 





CO(l-O) 



\ \ b \ — 

• • •• 

* • i 

■■* ••■••■*■••■•• ♦•-■■ 



CO(3-2) 




..• • .••. 



C4 


V 


CO(l-O) 


••••■•*-,• 


.•• 


•1*1 


..•• 


: ft 


^ 

....;..:. 




CO(3-2) 

-rV. 


• •• 

•♦^•.^, 

• • 



-400 -200 200 400 

velocity (km s M 



-400 -200 200 400 

velocity (km s M 



-400 -200 200 400 

velocity (km s M 



-400 -200 200 400 

velocity (km s M 



Figure D2. lit'-spectra of bin C galaxies. For each source, the top row (red) shows the CO( J = 1^0) hne and the bottom row (blue) 
shows the CO (J = 3—^2) line. For each line, for each source, the top panel shows the vector-averaged Real (solid) and Imaginary (dotted) 
amplitudes (Real,Imag in mJy beam"-*^) and the bottom panel shows phase (ph, degrees) of the calibrated uv data versus velocity (km s"-*^). 
With the higher resolution of the CO( J = 3^-2) data, we observe the peak of the CO( J = 3—^2) emission to be offset (< 1.5'' ~ 6 kpc) 
from the centers of galaxies CI and C2. In these cases, we calculate the uv spectra at the peak of the CO( J = 3—^2) emission. For the 
CO( J = 1 — )■ 0) line of galaxy CI, we average 2 channels together in order to increase the signal to noise ratio. 

coincident with noise- hke Imaginary amphtudes and phases of ~ 0°) over multiple velocity channels, indicative of a 
detection. We do not detect the CO(J = 3^2) line in galaxy CI. 

E. EXTRAORDINARY CO EMISSION IN GALAXIES A3 AND A8 

The integrated CO emission maps of galaxies A3 {z = 0.06) and A8 {z = 0.10) show evidence for disturbed 
morphologies. In both galaxies, we detect CO emission coincident with the optical emission, which shows a velocity 
gradient misaligned with respect to the optical major axis. In galaxy A8 we detect CO emission outside of the 
optical emission. Our observations are suggestive of CO emission outside the optical emission of A3 as well, but are 
not conclusive. In this Appendix, we investigate these emission components using iiv-spectra calculated at various 
positions in the field to distinguish true emission from noise. We discuss the more compelling case of galaxy A8 first, 
then A3. 

E.l. EGNoG A8 

In Figure El, we show the optical image, moment maps, integrated spectrum (as in Figure 3), and iiv-spectra of 
galaxy AS. In the top right panel, the solid blue (green dotted) line shows the integrated spectrum with (without) 
masking of the main CO component, indicated by the central white dotted ellipse in the moment map panels. The 
dashed red line shows the integrated spectrum of the emission component outside the optical galaxy (the offset white 
dotted ellipse in the map panels). The bottom three panels show the i^i;-spectra (amplitude and phase versus velocity) 
at 3 points. A, B and C, labeled on the map panels in the top row. 

Point A is centered on the CO emission coinciding with the optical emission. We only detect CO emission cor- 
responding to the upper right (north-west) half of the optical galaxy, at negative velocities (relative to the optical 
redshift from SDSS). The moment 1 map (intensity- weighted velocity) of this component shows a velocity gradient 
suggestive of a rotating disk. However, this velocity gradient does not lie along the optical major axis, but appears to 
be rotated by 45 to 90 degrees. 

Points B and C are at two positions along the long axis of the CO component that lies outside the optical galaxy 
(following the velocity gradient). The optical image from the SDSS shows no optical counterpart to this component 
of the CO emission. While the iiv-spectrum at point B appears dominated by noise, point C shows phases of ~ 0° 
over several velocity channels, indicative of real emission. This CO component shows a velocity gradient parallel to 
the major axis of the optical galaxy and is therefore not aligned with the velocity gradient of the CO disk component. 
However, the two components become roughly spatially coincident at negative velocities. 

In summary, the molecular gas in galaxy A8 is disturbed: CO emission is observed coinciding with the north-west 
half the optical galaxy but not in the south-east portion; the velocity gradient of the CO emission is not aligned with 
the optical major axis; and we observe significant CO emission outside the optical galaxy, with no optical counterpart. 
The optical image does not indicate interaction. We note that there is another galaxy in the field (seen on the left 
side of the optical image), but this galaxy is in the SDSS catalog with a photometric redshift of 0.28. The CO fiux of 
component outside the disk is 3.2 Jy km s~^, with ^center = —110 km s~^ and AV = 357 km s~^. This is very similar 
to the main CO component: ^co = 4.5 Jy km s~^, with ^center = —67 km s~^ and AV = 357 km s~^. If the emission 
outside of the optical disk traces dense molecular gas, the molecular gas mass of this component (assuming a Milky 
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Figure El. Optical image, moment maps, integrated spectra and ^t^'-spectra for EGNoG galaxy A8. The top four panels show (from left 
to right) the optical image, the moment map, the moment 1 map and the integrated spectrum (as in Figure 3) of the main part of the 
galaxy (solid blue shows spectrum with masking, dotted green without). The integrated spectrum of the external CO component is also 
shown by a red dashed line. The bottom three panels show the ^i^'-spectra (amplitude and phase) versus velocity of each of the 3 points 
(A,B and C) marked in the moment maps in the top row. 

Way-like conversion factor) would be 6.5 x 10^ Mq, roughly 70% the molecular mass of the component coincident with 
the optical emission. 

E.2. EGNoG A3 

In Figure E2, we show the optical image, moment maps, integrated spectrum (as in Figure 3), and ix'^-spectra of 
galaxy A3. The top row shows moment maps and integrated spectrum calculated in the velocity channels of the CO 
emission corresponding to the optical emission of the galaxy. In the top right panel, the solid blue (green dotted) 
line shows the integrated spectrum with (without) masking of the main CO component, indicated by the central 
white dotted ellipse in the moment map panels. The dashed red line shows the integrated spectrum of the emission 
component outside the optical galaxy (the offset white dotted ellipse in the map panels). The middle row is the same 
as the top row, but for the velocity channels of the negative velocity component of the emission lying outside of the 
optical galaxy. In the middle right panel, the solid blue line shows the integrated spectrum of the emission component 
outside the optical galaxy and the dashed red line shows the emission corresponding to the main CO component. The 
bottom three panels show the iiv-spectra (amplitude and phase versus velocity) at 3 points. A, B and C, labeled on 
the map panels in the top row. 

Point A is centered on the CO emission coinciding with the optical galaxy. The moment 1 map of this component 
shows a velocity gradient suggestive of a rotating disk, but the gradient is misaligned with respect to the optical major 
axis (rotated 45 to 90 degrees from the major axis, as in galaxy A8). 

Points B and C are at two positions in the emission component that lie outside the optical galaxy: point B marks 
an area that is bright at positive velocities (93 km s~^) and point C marks an area that is bright at negative velocities 
(—377 km s~^). The optical image from the SDSS shows no optical counterpart to either component of the CO 
emission. The ii'i;-spectra at points B and C (bottom row, middle and right panels) are suggestive of real emission but 
do not show a strong signal. The integrated flux from the positive (negative) velocity component is 2.6 Jy km s~^ (1.2 
Jy km s~^), with '^center = 93 (-377) km s~^ and AV = 138 (138) km s~^. If the emission is real and tracing dense 
molecular gas, the positive velocity component has a molecular gas mass of 1.5 x 10^ M© and the negative velocity 
component has 7.0 x 10^ M©. These masses are both small compared to the molecular gas mass of the main CO 
component, 1.4 x 10^^ M©. 

The misalignment (with respect to the optical major axis) of the velocity gradient of the main component of the CO 
emission is suggestive of a disturbed morphology. Our observations also suggest there is CO emission outside of the 
optical galaxy, but we do not detect this emission at a significant level. 
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Figure E2. Optical image, moment maps, integrated spectra and w-u-spectra for EGNoG galaxy A3. The top four panels show (from left 
to right) the optical image, the moment map, the moment 1 map and the integrated spectrum (as in Figure 3) of the main part of the 
galaxy. The integrated spectrum of the external CO component is shown by a red dashed line. The middle four panels are the same as the 
top four, but for the negative velocity component of the external CO emission. The red dashed line in the far right panel of the middle 
row shows the spectrum of the main part of the galaxy for reference. The bottom three panels show the ^i^'-spectra (amplitude and phase 
versus velocity) at each of the 3 points (A,B and C) marked in the moment maps in the top and middle rows. 



